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ABSTRACT 

Aims. We describe the design, operation, and first results of a photometric calibration project, called DICF (Direct Illumination 
Calibration Experiment), aiming at achieving precise instrumental calibration of optical telescopes. The heart of DICE is an illumi¬ 
nation device composed of 24 narrow-spectrum, high-intensity, light-emitting diodes (LED) chosen to cover the ultraviolet-to-near- 
infrared specUal range. It implements a point-like source placed at a finite distance from the telescope entrance pupil, yielding a flat 
field illumination that covers the entire field of view of the imager. The purpose of this system is to perform a lightweight routine mon¬ 
itoring of the imager passbands with a precision better than 5 per-mil on the relative passband normalisations and about 3 Angstroms 
on the filter cutolf positions. 

Methods. Prior to installation, the light source is calibrated on a spectrophotometric bench. As our fundamental metrology standard, 
we use a photodiode calibrated at the National Institute of Standards and Technology (NIST). The radiant intensity of each beam is 
mapped, and spectra are measured for each LED. All measurements are conducted at temperatures ranging from 0°C to 25°C in order 
to study the temperature dependence of the system. The photometric and spectroscopic measurements are combined into a model that 
predicts the spectral intensity of the source as a function of temperature. 

Results. We find that the calibration beams are stable at the lO^"* level - after taking the slight temperature dependence of the LED 
emission properties into account. We show that the spectral intensity of the source can be characterised with a precision of 3 Angstroms 
in wavelength, depending on how accurately we are able to calibrate the wavelength response of the mononochromator. In flux, we 
reach an accuracy of about 0.2 - 0.5% depending on how we understand the off-diagonal terms of the error budget affecting the 
calibration of the NIST photodiode. We describe how with a routine < 60-mn calibration program, the apparatus is able to constrain 
the imager passbands at the targeted precision levels. 

Key words. Instrumentation: miscellaneous - Cosmology: dark energy 


1. Introduction 

The measurement of the dark energy equation of state with type 
la supernovae (SNe la; see for recent examples; Suzuki et al. 
(2012); Betoule et al. (2014)) sets very strong constraints on 
the accuracy of the flux calibration of the imagers used to in¬ 
fer the SN la luminosity distances. Indeed, mapping the relative 
variation in SN distances with redshift boils down to compar¬ 
ing fluxes measured at both ends of the visible spectrum. A key 
concern is therefore the control of the imager throughput as a 
function of wavelength. In practice, this means, on the one hand, 
controlling the intercalibration of the fluxes, measured in differ¬ 
ent passbands, and on the other hand, measuring and monitoring 
the passband shapes, in particular, the position of the filter fronts. 

The best technique of passband intercalibration today relies 
on observations of stellar spectrophotometric standards, whose 
absolute spectral energy distribution (SED) is known with accu¬ 
racy. The most reliable sets of standards are probably those that 
were established to calibrate the instruments mounted on-board 
HST (CALSPEC, 2000). The CALSPEC flux scale is defined by 
the SED of three hot DA white dwarfs, as predicted by a NTLE 
model atmosphere code, itself tuned to the Balmer line profiles 
of each dwarf (Bohlin, 2007; Rauch et al., 2013). The calibra¬ 
tion of the HST STIS and NICMOS spectrographs obtained by 


observing these fundamental standards is then propagated to a 
larger network of redder and mostly fainter stars (Bohlin, 2010; 
Bohlin et al., 2014). 

This approach is robust, since it relies on more than one fun¬ 
damental calibrator. It is also sound, since each CALSPEC re¬ 
lease can be traced back to a well-identified series of HST obser¬ 
vations and to a particular version of a stellar atmosphere model. 
However, even DA white dwarfs are complex objects, and the 
modelling systematics that affect the CALSPEC flux scale are 
difficult to estimate with precision. 

Most recent supernova surveys have been anchored on the 
CALSPEC network. This was often implicit because people 
were using the SEDs of historical bright standards recalibrated 
with HST (Bohlin & Gilliland, 2004a,b) when interpreting their 
magnitudes as fluxes (see e.g. Amanullah et al., 2010). A few 
projects, such as SNLS and SDSS, tried to explicitly tie their cal¬ 
ibration to specific CALSPEC calibrators (Holtzman et al., 2008; 
Regnault et al., 2009). The calibration error budgets published 
by these authors were dominated by uncertainties introduced in 
the metrology chain that links the standard star observations to 
the science images. Recent work by Betoule et al. (2013) has 
shown, however, that it is now possible to transfer the CALSPEC 
flux scale with uncertainties comparable to the CALSPEC inter¬ 
nal uncertainties (ctj.,- ~ 3 - 5 mmag). 
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To improve on the current cosmological measurements, fu¬ 
ture dark energy surveys will need to intercalibrate their pass- 
bands at a few per-mil, which is below the current estimates of 
the CALSPEC uncertainty budget. Since modern, ground-based 
imagers display repeatabilities of 1-3 mmag (Montalto et ah, 
2007), they can in theory be calibrated to this precision level. 
However, it is not obvious yet that this ambitious goal is within 
reach of the stellar calibration techniques presented above. 

In this context, alternate calibration strategies have been pro¬ 
posed by several groups. Most have in common that they rely 
on laboratory flux standards instead of stellar models. Over 
the past two decades, flux metrology has indeed experienced 
dramatic improvements, moving from 1-2% source standards 
to detector standards accurate at the 0.1% level. In the visi¬ 
ble, the metrology standards are generally silicon photodiodes, 
whose calibration can be traced back to the Primary Optical 
Watt Radiometer (POWR, Houston & Rice, 2006), an electri¬ 
cal cryogenic substitution radiometer maintained by the National 
Institute of Standards and Technology (NIST) as the official im¬ 
plementation of the optical watt. To characterise the photodi¬ 
odes provided to its end users, NIST maintains a sophisticated 
metrology chain (Larason & Houston, 2008) involving interme¬ 
diate light sources, notably the SIRCUS laser facility (Brown 
et al., 2006, 2000) and the Spectral Comparator Facility (SCF). 
With a calibrated photodiode in hand and the goal of calibrating 
another light detector, the end user has no choice but to build an 
intermediate light source to transfer the NIST flux scale to his 
own instrument. As a consequence, most alternatives to stellar 
calibration implement variations around the generic metrology 
chain presented in Fig. 1. 

Among those, at least two groups have proposed to revisit 
the historical measurements of Hayes & Fatham (1975) using 
modern flux metrology techniques. The goal is to compare the 
CAFSPEC and NIST flux scales directly using bright standard 
stars. A first team is developing the ACCESS rocket-borne 40- 
cm telescope (Kaiser et al., 2010). The instrument is a spectro¬ 
graph, sensitive from 350-nm to 1.7qm and calibrated directly 
at the SIRCUS facility. The other team, comprising NIST sci¬ 
entists, is building a spectrophotometer at the focus of a 10-cm 
telescope, calibrated with an artificial source, which will target 
stars in the 0 < V < 5 magnitude range (e.g. McGraw et al., 
2012). Both projects will rely on HST observations to bridge the 
gap between their magnitude range and the CAFSPEC magni¬ 
tude range. They will permit a direct comparison between the 
CAFSPEC and NIST flux scales, which is an invaluable piece of 
information. This does not, however, relieve future surveys from 
monitoring their own instrument, in particular, from measuring 
their passbands in situ and checking for a possible evolution in 
the long term. 

Today, the required accuracy on the positioning of the pass- 
band cut-offs is as low as a fraction of a nanometre. As an ex¬ 
ample, for SNFS, decreasing the uncertainty on r-band position¬ 
ing from 3-nm to 2A decreases the uncertainty on w as much 
as adding 200 SNe la in the Hubble diagram. Passband models 
are usually built using pre-installation test-bench measurements 
of the imager optical components. This is not entirely satisfac¬ 
tory because Alters may evolve over time. For this reason, nearly 
every modern survey has plans to build and operate a dedicated 
calibration source, able to follow the shape and normalisation of 
the imager passbands in real time. A precursor in this domain 
is the apparatus described in Doi et al. (2010). It consists in a 
movable light source, designed to inject quasi-monochromatic 
light into the optics of the SDSS 2.5-m telescope, and measure 



Fig. 1. Generic metrology chain from POWR, the official imple¬ 
mentation of the optical watt, which is maintained at NIST, to 
the telescope imagers. 


its effective passbands in situ. This system permitted us to unveil 
significant variations in the blue part of the imager u band. 

Stubbs & Tonry (2006) push the concept one step further by 
describing how an in-situ narrow-band illumination system can 
be used to transfer the calibration of a NIST photodiode to an as¬ 
tronomical imager. Implementations of the procedure sketched 
in their paper are the systems developed for the Mosaic im¬ 
ager on the CTIO 4-m Blanco telescope (Stubbs et al., 2007) 
and for the PanSTARRS GigaPixel imager (Stubbs et al., 2010). 
Another notable example, built along similar design principles, 
is the DECal system (Marshall et al., 2013) now installed in 
the dome of the Blanco Telescope to monitor and calibrate the 
DECam imager (DePoy et al., 2008). All these systems consist 
in a large diffusive screen placed in front of the telescope pupil 
and illuminated with monochromatic light, generated either by a 
lamp coupled to monochromators or by a tunable laser. So far, 
these projects have mostly focussed on in-situ measurements of 
the passband shapes. Stubbs et al. (2010) was able to perform 
a comparison between a stellar and an instrumental calibration, 
reporting a qualitative agreement of ~ 5% between both. 

dice' was conceived by members of the SNFS collab¬ 
oration, building upon the lessons learned working with the 
MegaCam wide held imager. It shares most of the goals of 
the projects described above, but differs in several points of 
its design. The primary goal of the DICE project is to imple¬ 
ment the metrology chain sketched in Fig. 1, i.e. to transport 
the calibration carried by a NIST photodiode to the MegaCam 
imager and compare this calibration to CAFSPEC. This in¬ 
cludes performing a spectrophotometric calibration of the light 
source using the NIST photodiode as a metrology standard 
(Step O in Fig. 1), measuring the relative normalisations of the 
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MegaCam passbands from series of calibration exposures taken 
with DICE (Step @), and applying this calibration to observa¬ 
tions of CALSPEC standards (Step @). 

This paper is the first in a series of three papers that will 
describe Steps O, ©, and Eor now, we focus on Step O. 
We present the design and implementation of the DICE calibra¬ 
tion sources (§2). We then concentrate on the spectrophotomet- 
ric characterisation of the DICE sources in §3, §4, and §5, on 
a dedicated test bench equipped with a photodiode calibrated at 
NIST. Spectrophotometry is notoriously difficult, and we had to 
develop specific methods to remain immune to the measurement 
systematics. In §6, we venture into the area that will be covered 
in the next DICE papers: we sketch briefly how a DICE source 
can be used to calibrate a broadband imager, and we estimate 
the precision of the calibration that can be obtained from DICE 
observations by propagating the Step O error budget established 
in this paper. We conclude in §7. 

2. The DICE light source 

The design of the DICE light source has been described in 
Juramy (2006); Barrelet & Juramy (2008), and Juramy et al. 
(2008). Two demonstrators have been built and installed. A first 
apparatus (SnDICE) was designed in 2007 for the MegaCam 
wide held imager (Boulade et al., 2003) which equips the 3.6- 
m Canada-Erance-Hawaii Telescope (CEHT). In 2011-2012, a 
second light source (SkyDICE) was built for the Sky Mapper im¬ 
ager (Keller et al., 2007) with an improved design, building on 
the lessons learned with SnDICE. A third source was also built 
and kept at LPNHE for long-duration test-bench studies. 

A key requirement that guided our design effort is that a cal¬ 
ibration instrument must be as simple and easy to maintain as 
possible. Also, it must intrinsically be as stable as possible, with 
additional built-in redundancies that allow checking for long¬ 
term drifts of the intensity delivered into the telescope pupil. As 
a rule of thumb, we consider that to calibrate an instrument at 
the permil level, we must be able to measure and study any l%c 
drift of its response with a precision of 10% or better. This sets 
an ambitious stability goal on the light source of about 10 In 
the rest of this section, we present the most important aspects of 
the light source design. 

2.1. Calibration beam 

Ideally, a calibration light source should mimic the science ob¬ 
jects under study as much as possible. Since a supernova sur¬ 
vey is dealing primarily with point sources (supernovae and held 
stars), we should try to generate quasi-parallel beams, covering 
the entirety of the primary mirror (see figure 2(b)). Such a beam 
would result in a spot on the focal plane, and we could use the 
photometry code used in the survey photometry pipeline to es¬ 
timate its flux, thereby avoiding the systematic errors that arise 
from using different flux estimators. 

Unfortunately, building a good artificial star turns out to be 
difficult. We therefore deliberately opted for a different design 
(figure 2(a)). SnDICE is a point source located in the dome a few 
metres away from the telescope primary mirror and close to the 
object plane. The source generates a conical, quasi-Lambertian 
beam, of aperture ~ 2°, so slightly larger than the telescope’s 
angular acceptance. Such an illumination results in an almost 
uniform focal-plane illumination. 

As shown in figure 2, the calibration beam is radically dif¬ 
ferent from the science beam. In particular, the angular distri¬ 
bution of the light rays that hit the various optical surfaces (e.g. 



Mirror 


(a) DICE beam (b) Science beam 

Fig. 2. Left: the DICE calibration beam. A DICE light source 
consists in a quasi-point source, located at a finite distance from 
the telescope pupil. Such a beam generates a quasi-uniform illu¬ 
mination of the focal plane. Right: the science beam. The light 
from a distant source, such as a star, may be approximately mod¬ 
elled by a plane wave that Alls the entire telescope aperture. Such 
a beam is focussed on a single point of the focal plane. 


the interference Alters) is not comparable. However, this specific 
calibration beam has at least two very nice properties. Eirst, no 
intermediate elements (folding mirror, window, screen, etc.) are 
present between the light emission zone and the primary mirror. 
As a consequence, the design stays very simple, and the system 
is expected to be stable in the long term. Second, the structure 
of the beam is much simpler than the science beam, in the sense 
that each pixel sees photons that came through a unique path. In 
other words, there is a one-to-one relationship between the fo¬ 
cal plane elementary surface elements and the calibration beam 
elementary solid angles. It is therefore quite simple to predict 
the focal plane illumination, once one knows the beam radiant 
intensity map: the former follows from the latter from purely ge¬ 
ometrical considerations, involving propagation of light in free 
space, and through the optics. 

2.2. Light emitters 

Narrow-spectrum light emitting diodes (LED) were chosen as 
light emitters. As shown later in this paper, LEDs are extremely 
stable, as long as they are fed with stable currents and operated at 
a stable temperature. It it is relatively easy today to build current 
sources that are stable at a few 10“^ over a temperature range 
of a few degrees, and the LED emission properties vary with 
temperature in a smooth and simple manner. One of the purposes 
of this paper is to show that one can build a LED-based light 
source, delivering beams whose stability can be controlled at the 
level of a few 10 over long durations. 

LEDs do not emit monochromatic light. The typical EWHM 
of a LED spectrum is about 6A/A ~ 5 - 7% (i.e. 20-nm to 50- 
nm). This means that we need 20 to 25 LEDs to cover the entire 
visible spectrum - from 350-nm to 1100-nm. In figure 3, we 
show the sampling that could be obtained with the first proto¬ 
type, which was built for MegaCam (figure 3(a)). We also show, 
for comparison, what could be achieved four years later with 
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Fig. 3. Lower panels: (Left) coverage of the MegaCam filters with SnDICE (2008). The central wavelengths of all LEDs are dis¬ 
played with arrows. The shaded regions show the wavelength extension of each LED SED -within 7% of the spectrum peak. In 
2007, no narrow spectrum LED with sufficiently high power was available around 700-nm and above 900-nm. (Right) coverage 
of the SkyMapper filters with SkyDICE (2010). The position of about half of the LEDs were chosen so that they constrain the 
filter cutoffs well. The other LEDs were selected to fill the gaps between filter cutoffs with about one LED every 20/30-nm. Upper 
panels: (Right) sampling from the near-UV to the near-IR achievable with the LEDs listed in the OSRAM and Roithner catalogs. 
(Left) optimised sampling of the MegaCam passbands that can be implemented today with the LEDs available on the market. 


SkyDICE, our second prototype, built to calibrate SkyMapper 
(figure 3(b)). As can be seen, the diversity of LEDs available on 
the market improved very significantly in a few years. Today, by 
combining the catalogues of the three main LED manufacturers, 
it is theoretically possible to cover the entire spectral range of sil¬ 
icon imagers with about one LED every 10-20 nm (upper panel 
of figure 3(b)) or to design some sort of “optimal sampling” of 
the passbands (upper panel of figure 3(a)). 

In this design, we chose to sacrifice wavelength precision 
in favour of high-stability illumination. This makes sense, since 
what one actually needs is a follow-up more than an absolute 
measurement of the filter cutoff positions. The filter transmis¬ 
sions are measured well prior to installation, and our goal is in¬ 
stead to monitor any drift over the life time of the instrument. 
One of the aims of this paper is to estimate how precisely one 
can locate a filter front with such an instrument, which delivers a 
coarse wavelength sampling. We show in section 6 that excellent 
precision may be attained with a small number of exposures, as 
long as we can secure a precise model of the spectral intensity 
delivered by the light source. 

As temperature increases, the LED emission efficiency drops 
by about 0.5%/°C, and the mean wavelength of the emitted light 
shifts redwards by as much as 0.1 A/°C. We often refer to these 
variations as the “cooler-brighter” and “cooler-bluer” effects, re¬ 
spectively. As discussed in §3.2, these variations are generally 
linear and extremely reproducible. As a consequence, once each 
emitter has been characterised well, one only needs to imple¬ 
ment a real-time follow-up of the source temperature to account 
for these effects. 

The temperature variations induced by the LED itself are 
small, if detectable. Around 2007-2008, the typical power con¬ 
sumption of our LEDs operated at a low regime was of 50 to 
100 mW. With the new generation of LEDs available today, it is 
closer to 10 mW. A large fraction of this power is dissipated as 


heat, but it is easy to build a heat sink able to absorb all of it. 
We have verified that if a LED is correctly glued to a radiator 
with a heat-conducting glue, no noticeable temperature eleva¬ 
tion of the radiator itself can be detected. The only exceptions so 
far are less powerful LEDs (older models or UV LEDs), which 
have to be operated near their maximum power. Eor those emit¬ 
ters, we found typical radiator temperature elevations of about 
0.1 °C, with a maximum of 0.2°C for a UV LED, after an hour 
of operation). 

2.3. Light source design 

The mechanical design of the light source is illustrated in figure 
4. SnDICE and SkyDICE have very similar designs. SnDICE 
is a 150 mm x 150 mm x 300 mm modular box, made of eight 
almost identical anodised aluminium blocks, each pierced with 
25 apertures to let the light through. The LEDs are located 
on the back of the device, about 260-mm from the front face. 
The calibration beams exit through (2)9-mm apertures located 
on the front face of the device. This design permits conical 
2°wide beams to be generated. SkyDICE, the second prototype, 
is shorter (150 mmx 150 mmx227 mm) because it needs to gen¬ 
erate wider (3°) beams to cover the larger (5.7deg^) field of view 
of the imager. The modularity of the design allows us to adapt 
the aperture of the beam easily as a function of the imager field 
of view. 

The light source implements 24 calibration channels, each 
generating a conical Lambertian beam, in order to cover the 
350 nm < T < 1100 nm spectral range as evenly as possible. 
The LED currents are chosen so that a beam generates about 
1000 photoelectrons per second in each pixel. As the typical 
solid angle subtended by a pixel is about 10 sr, this means that 
the radiant intensity of the LED should generate of the order of 
10^^ y/s/sr, which corresponds to about 0.5 mW/sr at 500-nm. 


4 




























N. Regnault et al.: The DICE Calibration Project 



masks 


webcam 


LEDs 


radiator 


photodiodes 


front-face 


Fig. 4. Design of the DICE LED head. The LED head is made of 
8 aluminium blocks pierced with 24 holes corresponding to each 
of the 24 LED channels. Black masks are placed between each 
block to shape the beam and minimise the stray light. Each chan¬ 
nel is equipped with an off-axis control photodiode that monitors 
the light delivered by the LED in real time. The photodiodes are 
placed on a board just behind the front face of the device. The 
LEDs are in thermal contact with an aluminium radiator, in or¬ 
der to dissipate waste heat and monitor the LED temperature. 
An additional central 25th channel delivers a pencil beam used 
to control the relative orientation of the device with respect to 
the telescope. 



Fig. 5. Sketch of a calibration channel showing the emitting LED 
and the off-axis monitoring photodiode. 


The source is attached to the dome of the telescope. It can 
rotate around an altitude and an azimuth axis, in order to con¬ 
trol the alignment of the LED beams with the telescope axis. By 
moving the source and the telescope simultaneously, it is possi¬ 
ble to illuminate nearly every region of the primary mirror. 

A special LED channel, called the artificial planet is used to 
control the relative alignment of the device and the telescope. 
The planet light is generated by a wide spectrum LED cover¬ 
ing all the hlter passbands. The planet channel is equipped with 
a small ©10-mm convergent lens that transforms the isotropic 
LED beam into a quasi-parallel beam. Planet exposures produce 
a spot on the imager focal plane, along with ghosts, owing to in¬ 
ternal reflexions between the optical components. The position 
of the planet spot is a direct measurement of the angle between 
the telescope optical axis and the planet beam. It is measured 
with a precision of about 3 "(the EWHM of the spot itself being 
of about 20 "). This precision is slightly better than the precision 
of the motor axis encoders (0.0025 degree, 9 "). 

2.4. Generation of the LED currents 

The intensity of the light emitted by each LED is a function of 
the current injected into it. The LED currents are generated (and 
monitored in real time) by a custom-made backend board, con¬ 
nected to the illumination system with an analogue link. We at¬ 
tempted to build a current source that achieves an electrical sta- 



I 0.04 
"S 0.03 
g 0.02 

I 0.01 

3 0 

D-0.01 
-- 0.02 
■|-0.03 
£-0.04 


Fig. 7. Relative variations in the LED current, measured as the 
voltage drop across a resistor in series with the LED and mon¬ 
itored over a duration of 3 weeks. During this period, the LED 
head was operated without interruption on our spectrophotomet- 
ric test bench. During this run the temperature varied by about 
5°C. As can be seen, the current varies by a few 10“^ at most. 
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bility better than 10“^ over a temperature range of a few degrees. 
It is classically implemented as a transistor current source driven 
by a programmable voltage level and stabilised with a negative 
feedback loop (see e.g. Horowitz & Hill, 1989, p. 181). Eor re¬ 
dundancy, the feedback voltage level is also sampled at a rate 
of a few kHz by a ETC 1608 16-bit ADC and logged for offline 
checks. 

The accuracy and stability of this current source depends en¬ 
tirely on a few high-quality commercial components, used in the 
board, namely the voltage references of the ADC and DAC, the 
voltage division chains of the ADC and DAC, and the serial re¬ 
sistors. Extensive, long-duration tests of the generated currents 
show that they are stable at the 10^^-level (hgure 7). As the char¬ 
acteristics of the components may have a small temperature de¬ 
pendence, the temperature of the backend board is monitored in 
real time during data-taking. 

2.5. Redundancies 

Redundant checks are an essential part of the design. No mat¬ 
ter how stable the source is on a test bench, we need to prove 
that the light actually delivered during the calibration runs is not 
affected by unexpected fluctuations or long-term drifts. Eor this 
reason, several critical quantities are monitored in real time dur¬ 
ing operations. 

The temperature of the LEDs is the most important of all of 
these. What is measured in practice is the temperature of the ra¬ 
diator on which all LEDs are glued with a thermally conductive 
glue. Eor this purpose, we use a PTIOOO thermistor glued to the 
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radiator. What we monitor with this probe does not strictly cor¬ 
respond to the temperature of the LED junctions. It is, however, 
a surprisingly good proxy, which allows us to empirically stan¬ 
dardise each LED on a test bench and use these standardisation 
relations later during operations. 

Another critical point is the stability of the current source 
itself. We characterise it with two observables. Eirst, the tem¬ 
perature of the backend board, which may have an impact on 
the current generator, is monitored using a DS600 temperature 
sensor, mounted on the board itself Second, the current actually 
delivered to the LED (in practice, the voltage drop across a re¬ 
sistor mounted in series with the LED) is also logged for offline 
checks. 

Einally, we directly characterise the light actually delivered 
by the LED by placing an off-axis monitoring photodiode in 
each LED channel, close to the exit hole (see figure 5). We use 
5.8 mm x 5.8 mm Centronic^ OSD35-7 photodiodes, covering 
the full wavelength range of the illumination system with good 
efficiency in the UV. These control photodiodes are mounted on 
a board located just behind the front panel of the instrument (see 
figure 4). 

All these quantities are digitised on the backend board at 
a frequency that can be tuned (from IkHz to 32 kHz), so that 
one can study the possible sources of noise over a large band. 
The digital samples are stored in a 16MB buffer that is read out 
on demand by the DAQ system (typically after each calibration 
frame). 

Another miniaturised system, comprising a cooled, large 
area photodiode coupled to a ultra-low noise current amplifier 
has been built to monitor the light as it goes through the tele¬ 
scope optical path. It has not been used extensively in our anal¬ 
ysis. We describe it in Appendix E. 

2 . 6 . Data acquisition 

The DAQ and control system of these light sources run on a dedi¬ 
cated industrial PC (PC 104), which can be installed a few metres 
away from the source. It is connected to the backend board with a 
USB link and to the telescope control system (TCS) with an eth- 
ernet link. In its current implementation, it hosts a lightweight 
server, which interacts with the TCS through a variant of the 
XML-RPC protocol. The server communicates with the backend 
electronics and relays the orders sent by the TCS. In particular, 
it controls the LED current and the LED head orientation mo¬ 
tors, and it retrieves on demand the monitoring data stored on 
the backend board and sends them to the TCS so that they can 
be stored along with the calibration frames. 

2 . 7 . Operations 

All calibration frames are taken during daytime in order not to 
interfere with the telescope observing schedule. This requires 
the dome to be dark, or requires monitoring the ambiant lumi¬ 
nosity with the imager itself, complemented with an external de¬ 
vice, such as one of the modules described in Appendix E. In 
practice, we have found that in the years 2008-2010, the dome 
of CEHT was very dark with a contamination lower than 0.05 
ADU/s/pixel. Eor a 10 s exposure, with a typical level of ~ 5 000 
to 10 000 ADU, this yields a relative contamination by ambiant 
light of 5 10“^ to 10 at most. The enclosure of SkyMapper is 
slightly less light-tight, but the ambiant luminosity is also lower 
than 0.1 ADU/s, giving relative contaminations of a few 10“"^. 

^ http://www.centronic.co.uk 


With such low contamination levels, we did not implement any 
correction at the pixel level. During a later series of calibration 
runs performed at the CEHT from January to August 2014, the 
dome was found to be significantly brighter, following the in¬ 
stallation of venting apertures (several ADU/s/pixel). The am¬ 
bient luminosity was then monitored by interlacing the calibra¬ 
tion frames with dark dome exposures taken in the same position 
with the LED turned off. This nearly doubled the duration of a 
typical calibration run from ~ 50’ to 1 h 40’ (see Table 2 of §6.3). 

A calibration run goes as follows. The telescope and dome 
motion are decoupled, and the telescope points inside the dome 
in the direction of the illumination device. The optical axes of 
both instruments are then aligned using the planet beam (see 
section 2), as a guide. Once the alignment of both instruments 
is known, within an acceptable range, series of calibration expo¬ 
sures may be taken using the main LEDs. 

Several types of calibration exposures may be taken. The sta¬ 
bility of the readout electronics is studied with repeated expo¬ 
sures of the same LED, taking advantage of the 10 stability of 
the illumination system. The linearity of the imager is checked 
using illumination ramps (i.e. exposures of the same LED of 
longer and longer exposure time). Einally, one measures the in¬ 
strument passbands with series of calibration frames taken with 
all the LEDs matching the passband under study. 

As we see in section 6, calibrating one single filter requires 
taking about four to eight exposures (each with a different LED), 
of one to ten seconds each. The length of a calibration sequence 
is therefore dominated by the imager readout time (40 s for 
MegaCam) and never exceeds ten minutes. As a result, it is 
possible to check the five to six filters that equip MegaCam or 
SkyMapper in less than one hour. 

3. Spectrophotometric test bench: Overview 

3.1. Goai 

The primary goal of the test bench studies is to transfer the cal¬ 
ibration carried by the NIST photodiode to the light source (ar¬ 
row labelled O in figure 1). In practice this means that we want 
to characterise the emissivity of each LED using the NIST pho¬ 
todiode as our primary standard. 

The physical quantity which describes the emissivity of a 
point source is its spectral intensity, S{A,u). It is defined as the 
power emitted by the source, per unit wavelength and per unit 
solid angle (at a given wavelength A and in a given direction 
u). Since the properties of many types of LEDs are sensitive to 
temperature, we expect S to also be a function of T: >S(/l, T, u). 

Our goal is to build a smooth model S (T, T, u) of the LED’s 
true spectral intensity >S(/1, T, u), which is valid in the tempera¬ 
ture range 0°C < T < -i-25°C and which is typical of what is 
measured in most telescope enclosures. Along with this model, 
we would like to build an uncertainty model to account for the 
finite precision of the test-bench measurements and the finite sta¬ 
bility of the source. These two ingredients will then allow us to 
predict, with a known accuracy, the amount of calibration light 
delivered by the source for any real calibration exposure. 

To build such a model, we need to accumulate enough mea¬ 
surements at temperatures typical of what is measured in a tele¬ 
scope enclosure. We also need to characterise the short-term and 
long-term stability of the source. In the next sections, we de¬ 
scribe the test-bench configurations, the measurements that are 
taken, and the way we combine them to build a smooth estimate 
of the LED spectral intensities. 
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Fig. 8. Spectro-photometric test bench. Left panel: photometric bench. The photodiode is placed on a linear table and moved in 
planes parallel to the z-axis to map each LED beam. Right panel: a monochromator is inserted between the source and the photodi¬ 
ode. 


3.2. Test-bench setup 

A sketch of the calibration test bench is shown in hgure 8. The 
light source is placed on one end of a 2.5-m long optical table 
and mounted on an X-Y Kinetic Systems KVP-lOO linear table, 
with a total stroke length of 150-mm. 

A 1 cm^ Hamamatsu S2281 Si photodiode was purchased 
from NIST to be used as our primary standard. Following the 
NIST procedures, it is operated unbiased at room temperature^. 
Depending on the bench setup (see below), the photocurrents 
generated by the photodiode vary from a few nano-amperes (in 
photometric mode) down to a few pico-amperes (in spectro¬ 
scopic mode). They are measured with a Keithley 6514 feedback 
picoammeter. 

The calibrated photodiode is mounted on a KVP-lOO X-Y 
linear table that allows for a stroke length of about 300-mm. A 
third motorisation permits moving it by about 1500-mm along 
the Z-axis. The repeatability of the KVP-lOO tables is of 1 ^im, 
allowing us to control the relative positioning of the source and 
the detector very precisely. 

The test bench is placed ina2mx2mx3.5m dark enclo¬ 
sure. The enclosure walls are insulated, and using a powerful air 
conditioning system, we manage to cool it down to a tempera¬ 
ture of ~ 0°C. The bench is not strictly speaking thermalised, 
since its temperature is neither regulated nor perfectly uniform 
within the whole volume. We compensate for this by monitoring 
the temperature of the key parts of the test bench using PTIOOO 
thermistances or type-K thermocouples. 

3.3. Assumptions 

S(A, T, u) depends on four scalar parameters. This means that 
to pave the full parameter space, we need to accumulate a large 
number of spectra. However, we have verihed on six LED mod¬ 
els, chosen so as to cover the full spectral range and the full range 
of available LED technologies, that the spectra are essentially in¬ 
dependent of the direction of emission u (see Appendix D for a 
full discussion). As a consequence, we assume that the spectral 
intensity of all the LEDs that equip the DICE light sources can 
be written as 

S{A,T)xB(u) (1) 

where S is the spectral intensity of the LED in a (arbitrary) ref¬ 
erence direction, while S(h) is a dimensionless function, which 
accounts for the variations in the beam intensity as a function 
of the angle of emission. Here, tB{u) is normalised to one in the 
reference direction. 

^ This is sometimes called “photovoltaic mode” as opposed to the 
photoconductive mode where the photodiode is reverse-biased. 


In what follows, we refer to .S(/l, T) as the LED “spectrum”, 
keeping in mind that it is actually a spectral intensity. The di¬ 
mensionless quantity S(h) is called the “beam map”. 

3.4. Photometric and spectroscopic measurements 

The beam maps S can be measured simply by intercepting the 
beam with a calibrated photodiode placed at a known distance 
from the source and moved with respect to the source, in or¬ 
der to sample the whole beam. The measurement of is a little 
more complex. We need to perform spectroscopic measurements 
by inserting a monochromator between the source and the cali¬ 
brated photodiode. 

The calibration of the light source is therefore performed in 
two distinct steps. First, we simply map the radiant intensity of 
each calibration beam by moving the standard photodiode in a 
series of planes orthogonal to the Z-axis. These calibration se¬ 
quences, called hereafter “photometric calibration sequences”, 
are performed at about 10 to 15 different temperatures, ranging 
between ~ 0°C and ~ 25 °C. They allow us to study how the in¬ 
tensity delivered by each LED varies with temperature and how 
the intensity varies with the direction of emission (beam map). 
They also permit us to assess the stability of the source. The 
analysis of this dataset is discussed in section 4. 

In a second step, we insert a Digikrom DK240 Czemy- 
Turner monochromator between the light source and the cali¬ 
brated photodiode. The LEDs are positioned in turn in front of 
the monochromator entrance slit, while the photodiode senses 
the intensity that comes out of the exit slit. These measure¬ 
ments are performed at about ten distinct temperatures between 
0°C and room temperature. This spectroscopic dataset, com¬ 
bined with the photometric measurements described above, al¬ 
lows us to derive smooth models of the LED spectral intensity. 
The spectroscopic measurements are described in section 5 with 
additional details given in appendix C. 

4. Photometric calibration of the light source 

The goal of the photometric calibration studies is to deter¬ 
mine the absolute normalisation of the LED spectral intensities 
>S(/1, T) and also to measure the beam maps S(h) (see section 
3.3). In this conhguration of the bench, the calibrated photodi¬ 
ode samples directly the beam light (left panel of hgure 8). The 
photocurrent registered with the Keithley 6514 picoammeter can 
then be written as the product of two simple quantities: 

^Iphot “ !B(u)xJ(T) (2) 

where !B{u) is the beam map dehned in section 3.3, T is the 
LED temperature (or at least a proxy for it), and J'iT) the pho- 
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tocurrent generated when the photodiode is placed at a specific 
reference position ro with respect to the source. At this reference 
position, S is conventionally set to 1, the photodiode subtends a 
solid angle dflo, and we have 

J(T) = dQo X J 77(T) *S(T, T) dA, (3) 

where //(T) is the photodiode efficiency reported by NIST. The 
main difficulty is to control the relative positions and orientations 
of the photodiode with respect to the source. 

Since acquiring detailed beam maps is time consuming (~ 30 
to 60 minutes per map), we optimise the photometric measure¬ 
ments as follows. In a first series of measurements, we concen¬ 
trate on a few specific beam locations, keeping the photodiode 
fixed, while varying the bench temperature. These sequences, 
called “minimaps” are much faster to acquire (about 260 sec¬ 
onds) and yield about ten independent measurements of each 
selected beam location (in particular, the central region of the 
beam, taken as a reference). They are targeted at measuring 
J^{T) for each LED, i.e. the relative variations in the LED emis¬ 
sion with temperature. The analysis of the minimaps is presented 
in section 4.1 below. 

We also realise fine-grained maps of the calibration beams at 
two or three temperatures, keeping the temperature of the bench 
as constant as possible during data taking. Erom this data, we 
obtain the !B(u) maps described above, and we verify that these 
maps are stable with temperature (see section 4.3). Such mea¬ 
surements are taken at several distances to the source, in order to 
verify the projectivity of the beam against scattered light and to 
check that we control the LED-source geometry well. 


LED D1 



Fig. 9. Top panel; Cooler-brighter relation measured from 13 
minimaps for a Roithner-Lasertechnik APG2C1-690 emitter, 
mounted on SkyDICE. Bottom panel; residuals. The small black 
error bars display the uncertainties on the minimap measure¬ 
ments (combined). The larger grey error bars take the measure¬ 
ment repeatability between maps into account (i.e. over time 
scales of one hour). 


4.1. Minimaps: LED emission versus temperature 

In most cases, the LED emissivity decreases with temperature. 
We refer to this as the “cooler-brighter” relation already men¬ 
tioned above. It generally obeys a linear law; 


"cooler-brighter'' slopes (SkyDI CE) 
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J'{T) - obc X (Tled - To) -H J'o (4) 

where J^(T) is the photodiode current defined above,?’ a proxy 
for the LED temperature, typically the temperature of the LED 
radiator, and To is an arbitrary temperature pivot. 

Ten to fifteen minimaps are taken for each LED at tempera¬ 
tures ranging from 0°C to 25°C. The cooler-brighter law is fitted 
on the measurements of the central region of the beam. An exam¬ 
ple of the fit, along with the residuals to the linear law, is shown 
in figure 9 for one of the LEDs that equips SkyDICE. 

Our determinations of the LED radiant intensities at 25°C 
and of the cooler-brighter slopes are reported in Tables A. 1 and 
A.2 in the appendix for SnDICE and SkyDICE, respectively. 
Whenever possible, the SkyDICE LEDs have been tuned to de¬ 
liver of the order of 0.1 mW/sr. On the test bench, placing the 
1 cm^ NIST photodiode about 2 m away from the source, this 
translates into typical photodiode currents of a few 10^ pA, as 
shown in figure 9. Such a current level is easily measured with 
a picoammeter. In figure 10 we summarise the intensity of the 
cooler-brighter effect for all the LEDs mounted on SkyDICE. 
On average, it is slightly lower than 0.5%/°C for all LEDs (re¬ 
gardless of the LED technology), except for two blue InGaN 
emitters that can reach over 1%/°C. We notice that two LEDs 
emitting around 500-nm actually exhibit a warmer-brighter re¬ 
lation. It may be because the nominal currents that were chosen 
for these LEDs are below the normal range recommended by the 
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Fig. 10. Cooler-brighter slopes for the SkyDICE LEDs. The var¬ 
ious LED technologies are coded with different symbols (see 
Appendix A for details). The amplitude of the effect is similar for 
all LEDs at about 0.5%/°C. Eor two LEDs, which are operated 
at very low currents, we actually observe a “brighter-warmer” 
relation. 


vendor. What matters for our application is that the effect is well 
measured and reproducible. 

Only one of the LEDs mounted on SkyDICE clearly displays 
small deviations from a linear cooler-brighter law. It is a Golden 
Dragon® LED of type LD W5AM, emitting at ~ 450-nm. The 
intensity of the effect is small (< 0.2% peak-to-peak over the full 
temperature range), and we model it with a second-order poly- 
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Fig. 11. Upper panel: variations in the source temperature over a 
3-week run designed to test the system stability. Middle panel: 
black triangles: relative variations of the raw flux delivered by 
the source. Blue lower triangles: relative variations in the flux 
after correction for the cooler-brighter effect. Red upper trian¬ 
gles: relative variations in the flux after correction by the control 
photodiode measurements. Lower panel: zoom on the corrected 
flux baselines. 


nomial. For all the other LEDs, it is not possible with the data 
in hand to distinguish between small deviations from linearity 
from the bench or source instabilities. 

Special care has been taken into evaluating the measurement 
repeatability empirically on various time scales. The black error 
bars in the lower panel of figure 9 display the empirical variabil¬ 
ity of the central beam region, measured on each single minimap 
(i.e. over about 260 seconds). It is about 0.1 pA (about 5 nW/sr) 
and depends slightly on whether the cooling system is on or off. 
When combining the ten minimap measurements together, it rep¬ 
resents a negligible contribution to the flux uncertainty: about 
1.5 nW/sr, for a nominal flux of about 0.1 mW/sr, i.e. a few 10“^. 
The grey error bars in the same figure represent the map-to-map 
variability. This contribution depends on the LEDs, so we at¬ 
tribute it to some longer term-variability of the source, and to 
a minor extent, to some additional variability of the bench. We 
model it with a noise pedestal (Tint, adjusted iteratively to obtain 
a reduced of unity. Eor most LEDs, it represents a little less 
than 5 10'"^ of the LED nominal flux. Eor the two faintest UV 
LEDs, it is slightly larger (2 10“^). 

4.2. Long-term stability of the DICE illumination device 

In the previous section, we have shown, from the study of the 
residuals to the brighter-cooler law, that the source is stable at the 
level of a few 10 (rms) over the duration of a calibration run 
(i.e. 24 hours). Since the stability of the light source is a critical 
aspect of the design, it has been checked over longer durations 
(i.e. weeks). Eor these tests, the bench was configured in photo¬ 
metric mode, and repeated minimaps of each calibration chan¬ 
nel were acquired. The system cycles through all LEDs over and 
over again for a little more than three weeks. The tests were per¬ 
formed at room temperature, each LED being fed with its nomi¬ 


nal current. All the auxiliary quantities (LED temperature, back¬ 
end board temperature, control photodiode current, etc.) were 
logged, along with the NIST photodiode measurements. 

Eigure 11 shows the flux baseline that was acquired for a 
typical LED, along with the LED temperatures. Significant, anti- 
correlated variations of the temperature (~ 5°C) and the LED 
flux (~ 2%) are noticeable over the course of the run. This is the 
consequence of the cooler-brighter relation studied in the previ¬ 
ous section. 

We can take advantage of the cooler-brighter effect to predict 
the variations in the LED flux, since we monitor the temperature 
of the LED radiator Tled- We also measure the temperature of 
the backend board that generates the LED currents, Tbb, which 
has been found to correlate slightly with the LED flux for some 
channels. We use these two variables to predict the variations in 
the LED flux over the long run: 

^fpLED - flLED X ^Tled + ^bb X d^bb- (5) 

We find that Fled and Tbb permit parametrizing most of the vari¬ 
ations in the LED flux over time as shown in figure 11 . We show 
the residual dispersions attained over three weeks of continu¬ 
ous monitoring in figure 12. It is in the range 2.10“^ < < 

10“^ (rms) with a median value of 6 10 (rms). Only two LEDs 
are less stable than 5 10“^. One is a UV LED (340-nm) with a 
large threshhold voltage (about 5V). Such voltages are slightly 
higher than what our current source was designed for, and it 
seems to become less stable in this regime. Eortunately, the mod¬ 
ern UV LEDs on the market require lower threshhold voltages, 
so we should not have to deal with this problem in the future. The 
other unstable LED, a red APG2C1-760 distributed by Roithner- 
Lasertechnik, has displayed a true instability during the run, its 
flux jumping by almost 2% about two weeks after the beginning 
of the operations. This behaviour has not been explained yet and 
has not been seen with any other LED. It has no counterpart on 
the monitoring data, besides the LED control photodiode placed 
close to the exit hole (see Appendix B). 

We have developed an alternate and redundant standardisa¬ 
tion method that relies on the direct measurements of the LED 
fluxes performed with the off-axis control photodiodes. Since the 
photodiode currents are digitised on the backend board, we also 
include the backend temperature in the standardisation relation: 

^Vled - n/X 5(ppd -H x d^bb- (6) 

Again, we observe very low residual dispersions of the standard¬ 
ised flux (2 10“^ < < 10“^) (rms) with a median value 

of 4 10“^ (rms), as shown in figure 12. Again, there are a few 
outliers: the two LEDs discussed above, as well as a third LED 
whose control photodiode seems to be malfunctioning. These 
channels still display repeatabilities that are better than 1% over 
weeks. 

We plan to extend the duration of these intensive stability 
tests to longer durations of the order of a few months. In a near 
future, we will also re-calibrate the sources that are currently in 
operation (SnDICE at CEHT and SkyDICE at SSO). This will 
allow us to check the stability of their calibration over durations 
of a few years. 

As a conclusion, we have built a calibrated source that is in¬ 
trinsically stable at the 10^^-level. The fluctuations of the light 
delivered by each LED can be parametrised as a function of sim¬ 
ple auxiliary observables, notably the temperature of the LED 
radiator, Teed- Using two independent methods (one based on 
monitoring Teed and another based on the measurements per¬ 
formed with the off-axis control photodiodes), we have shown 
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Fig. 12. Red histogram: stability of the LED flux (i.e. the NIST 
photodiode current) after correction for the cooler-brighter ef¬ 
fect (and for the backend board temperature variations). Blue 
histogram: stability of the LED channels after correction for 
the control photodiode measurements. The two outliers above 
5.10“^ are LEDs in the blue histogram that correspond to chan¬ 
nels that are intrinsically unstable. The outliers in the red his¬ 
togram correspond to the same LEDs, plus one LED whose con¬ 
trol photodiode is faulty. 


that we can correct for the variations of the source at the level of 
a few 10 During operations, we generally choose to use the 
temperature method as our primary standardisation technique 
and to use the control photodiode for redundancy checks. The 
photodiode and temperature probes are read before and after 
each calibration exposure, therefore either of the two techniques 
can be used for each exposure. Since both methods give similar 
results, they can be interchanged in the future. 

Because of its stability, our light source qualifies as a “stan¬ 
dard light source”. Once calibrated, it may be used to dissemi¬ 
nate a flux calibration to remote locations. 

4.3. Beam maps: LED emission versus direction 

We now describe how we derive fine-grained beam maps (S) 
from the detailed photometric scans realised at nearly constant 
temperature (figure 13). We have found that the relative vari¬ 
ations in the LED emissivity with direction are small (1% at 
most) and smooth. They sometimes display a complicated shape, 
which does not suggest an analytic parametrisation. We therefore 
develop them on the basis of splines: 

S(u) = Y^PpBp(u): (7) 

p 

and we fit for the (3 parameters. The main source of systematics 
here is the determination of dQ, i.e. the control of the test-bench 
geometry. 

The stability of the beam maps themselves has been stud¬ 
ied in detail over long durations. The temperature fluctuations 
may indeed dilate the metal, hence slightly change the relative 
positioning of the LED with respect to the masks and pieces of 
baffles that shape the beam. We have found that the fluctuations 
of the map profiles are always lower than 5 10 


3D row-map of LED05 



Fig. 13. Beam map of LED LVW5AM (OSRAM, Golden 
Dragon® series) mounted on the SkyDICE device. 


5. Spectroscopic calibration of the light source 

We now turn to the determination of the LED spectra. This re¬ 
quires additional measurements, which are performed on our 
spectroscopic test bench (see right panel of figure 8). This bench 
is similar to the photometric test bench described above, except 
for the presence of a Czerny-Tumer type monochromator, in¬ 
serted between the LED source and the calibrated photodiode. 
In this configuration, the photocurrent generated by the photodi¬ 
ode is then given by 

/|spec = [r](4) ■ TmiA) ■ S(A, T)] ® WJA) (8) 

where rjiA) is the photodiode efficiency reported by NIST, T^iA) 
and WmiA) are the transmission and the spectral response of the 
monochromator, respectively. 

The key point here is the control of the monochromator cal¬ 
ibration. Erom the equation above, we see that we need to check 
(1) the wavelength calibration of the monochromator, (2) its 
transmission r„,(/l), and (3) its spectral response, Wm(A). This 
work is described in detail in Appendix C and briefly sum¬ 
marised in section 5.3.1 below. 

5.1. Modeiiing the LED spectra 

Our goal is to build a smooth model of the LED spectral intensi¬ 
ties as a function of wavelength and temperature. We choose to 
develop this model on the basis of two-dimensional B-splines: 

S(A,T)^YjepBp(A,T) (9) 

p 

where the Bp{A, T) functions are two-dimensional splines of or¬ 
der 3 to describe the wavelength variations, and order 2 to model 
the temperature variations. We use three nodes placed every 
10°C to describe the temperature and 0(100) nodes, placed ev¬ 
ery 2-nm to capture the spectral shapes. 

Because the uncertainty on the absolute normalisation of the 
monochromator transmission is difficult to assess, we decided to 
fit for the (unknown) normalisation fs of each spectrum measure¬ 
ment s. This means that the only piece of information we extract 
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Fig. 14. Upper panel; Comparison between a typical LED spec- Fig. 16. Cooler-bluer slopes measured on the SkyDICE spectra, 
tram (SkyDICE Golden Dragon® LV W5AM) and the htted 
model (red line). Bottom panel: residuals. 


from the spectroscopic measurements is related to the spectrum 
shapes and not the spectrum normalisations. Our model of the 
photodiode current (equation 8 above) becomes 

f I spec ~ fs ^ ^ Op [ 77 (d) ■ W) ■ BpiA, r)] ® w^iA). (10) 

p 

By construction, there is a perfect degeneracy between the ffpS 
and the f'^ s. We break it by requiring that the absolute normalisa¬ 
tion of the S models is determined by the photometric measure¬ 
ments (performed without a monochromator). In practice, this is 
done by adding the following term to the ht;^^^: 

J 77 (d) ■ 5(^,7,) 7/dj (11) 

where JL(L) is the photocurrent measured in a reference posi¬ 
tion, as dehned in section 4. With this approach, the uncertainty 
affecting the monochromator transmission is entirely absorbed, 
and the absolute normalisation of the LED spectral intensity 
models is set exclusively by the photometric measurements. 

We estimate the LED spectral intensities by htting the model 
above on the photodiode current measurements and marginaliz¬ 
ing on the fs nuisance parameters. 

5.2. Results 

The spectral intensity of the SnDICE and SkyDICE LEDs, >S(d), 
have been reconstructed using the method described above, com¬ 
bining the spectroscopic and photometric data. On average, 
about 15 spectra and a similar number of minimaps have been 
taken for each LED at temperatures ranging from 2°C to 25°C. 
In hgure 14 we show a comparison between the model and a 
typical spectrum. Eigure 15 shows the spectral intensity surfaces 
obtained for typical blue and red LEDs. These models sum¬ 
marise the behaviour of the source. Their normalisation is set 
by the photometric dataset. By construction they incorporate the 
cooler-brighter and cooler-bluer effects. Once we know the op¬ 
erating temperature and the distance between the source and the 
telescope aperture, we can predict the spectral intensity delivered 
on the primary mirror. 

As can be seen in figure 15, the cooler-bluer effect does not 
have the same intensity for all LEDs. In figure 16 we report the 


variations in the mean spectrum wavelength as a function of tem¬ 
perature (estimated from the model, for all the LEDs that equip 
SkyDICE). On average, the effect is about 1A/°C. A trend is 
clearly visible, indicating that the red LEDs are more sensitive 
to the effect than the blue LEDs. The two LEDs that have been 
found to exhibit a warmer-brighter behaviour (section 4.1) also 
follow a warmer-bluer relation. 

The spectral intensity models are affected by uncertainties 
owing to the hnite number of spectra and flux measurements. 
These uncertainties are statistical in nature, but behave as sys- 
tematics in any subsequent analysis relying on these spectral in¬ 
tensity models. Eor this reason, we are careful to propagate them 
in an exact way (i.e. including the off-diagonal terms). They are 
derived from the covariance matrix of the spectral intensity fit 
presented in the previous section. We hnd that with the datasets 
in hand, the relative uncertainties on the LED broadband fluxes 
predicted with the spectral intensity model are below 0.05%. 

5.3. Test-bench systematics 

We now examine the systematic uncertainties that have an im¬ 
pact on our spectral intensity models, S (A, T). We identify three 
main sources of systematics: the wavelength calibration of the 
monochromator, the transmission of the same monochromator, 
and the accuracy of the NIST photodiode efficiency t]{A). 

5.3.1. Monochromator 

Wavelength calibration The wavelength calibration of the de¬ 
vice is checked (at several temperatures), using a series of cali¬ 
bration lamps, notably a sodium lamp and a polymetallic lamp. 
Erom this dataset we obtain a correction to the calibration given 
by the manufacturer. This correction does not exceed 0.1-nm 
in amplitude. We hnd a small linear dependence of the wave¬ 
length calibration with temperature (see Table C.l) of about 
0.025 nm/°C, which is taken into account. The dominant system¬ 
atics affecting the wavelength calibration is related to the posi¬ 
tioning uncertainty of the source with respect to the monochro¬ 
mator entrance slit. We found that this contribution never ex¬ 
ceeds 0.1-nm. Given the limited precision of the alignment pro¬ 
cedure and the small number of calibration runs, we conserva¬ 
tively set the wavelength calibration systematics to 0.3-nm (see 
Table 1). 
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Fig. 15. Spectral intensity surfaces S (A, T) for two SkyDICE LEDs. The model incorporates the cooler-brighter and cooler-bluer 
effects. 


Table 1. Summary of the test bench systematics. 



Uncertainty 

(I 0 -) 

Comment 

Monochromator 


Wavelength calibration 

^Ebert 

fbiazc (grating #1) 
fbiaze (grating #2) 
fbiaze (grating #3) 

0.3 nm 

1“ 

0.17° 

0.28° 

0.43° 

{Aled) < 450 nm 

450 nm < {Aled) < 750 nm 
{Aled) > 750 nm 

Hamamatsu S2281 calibrated @ NIST 



From the calibration 

On 1ST 

2.9 10-5 nm-‘ 

uncertainties 

Pnist 

0.002 

provided by NIST. 


Transmission The transmission of the monochromator, T,„(A), 
is measured at several discrete wavelengths. Eor Czerny-Turner 
designs, it is relatively easy to compute the global shape of 
T„j(A), as a function of two specific angles: the so-called “Ebert 
angle” oe, which depends on the optical design of the device, 
and the “blaze angle” ebiaze|i of each grating i being used. We 
use these continuous models (one for each grating) fitted on 
the discrete transmission measurements as our estimates of the 
monochromator transmission (see figure C.l). At first order, 
T„j(A) only depends on the etiaze parameter. The uncertainty on 
the monochromator transmission is estimated by propagating the 
uncertainties on our estimates of the fbiaze angles. They are re¬ 
ported in Table 1. 

The method used to extract the LED spectrophotometric 
models (equations 10 and 11) ensures that the model normal¬ 
isation only depends on the photometric measurements (with¬ 
out the monochromator). This means that the uncertainties on 
the monochromator transmission have no impact on the LED- 
to-LED calibration and may only affect the spectrum shapes. 
To fix ideas, we find that varying the Ebert angle by one de¬ 
gree has an impact on the spectrum mean wavelength ((A) - 


f AS (A)dAf f S (A)dA), which is lower than 1 A. We also find that 
altering the monochromator transmission model around 550-nm, 
as suggested by the residuals of figure C.l, results in a shift in 
wavelength that is also lower than 1 A. 

5.3.2. NIST 

The calibration of the NIST photodiode is itself uncertain. NIST 
provides its clients with an error budget and reports uncertainties 
of about ~ 0.2% between 400-nm and 950-nm and up to 1% in 
the UV and near-IR (figure 17). We expect a fraction of the un¬ 
certainties affecting measurements at different wavelengths to be 
correlated. And since we are primarily interested in the relative 
calibration of our passbands, it is essential for us to account for 
these off-diagonal terms. As of today, we have not been able to 
obtain this information from NIST. We therefore built two differ¬ 
ent error models, depending on how we choose to interpret the 
NIST uncertainties. 

In the “best case scenario”, we assume that what is uncertain 
in the photodiode efficiency is primarily a grey scale. This means 
that the NIST uncertainties are positively fully correlated. Since 
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Fig. 17. Upper panel; efficiency (in A/W) of the Hamamatsu 
S2281 photodiode calibrated at NIST, purchased for this project. 
Bottom panel: black points: uncertainties reported by NIST (in 
percentage); dashed red line: diagonal elements of the “worst 
case” interpretation of the NIST error budget, as discussed in 
section 5.3.2; dashed blue line; diagonal elements of the “best 
case” interpretation of the NIST ERROR budget discussed in 
the same section. 


our goal is to measure the relative normalisation of the imager 
passbands, this would be an ideal situation, because the uncer¬ 
tainty on the global flux scale cancels out when comparing two 
different passbands. 

In the “worst case scenario”, we assume that there is a 
“colour uncertainty” affecting the efficiency reported by NIST. 
This can be modelled with two random variables aNisx and 
y^NiST, of variance ctq, and cr^, and such that 

?7nist(T) = 7/true(d) X (ctnist (d - T) h- ^nist + l) > (12) 

where cr„ and ojg are chosen so that the uncertainties com¬ 
puted from the equation above stay compatible with the error 
budget reported by NIST, and d is a pivot value (chosen to be 
d = 650 nm). Their values are listed in Table 1. 

This model is of course a little extreme because it assumes a 
correlation length that spans the entire visible range. This would 
mean that the uncertainty reported by NIST lies with POWR, the 
primary flux standard. This is not the case, since for the first steps 
of its metrology chain, NIST reports uncertainties of a few 10 
about one order of magnitude lower. It is more likely that the 
main fraction of the uncertainty budget is introduced in the last 
step at the SCF facility, and it is difficult to see how SCF could 
introduce error with an infinite correlation length in wavelength. 
However, in the present state of our knowledge, exploring these 
two extreme cases is the best we can do, because they bracket 
the true uncertainties. 

5.3.3. Propagating the bench systematics 

The test bench systematics affect how we reconstruct the LED 
spectral intensities {S ). They may, for example, shift the S (A) 
models in wavelength or distort their shape, and this must be ac¬ 
counted for. Regardless of the propagation method we choose, 
we need to quantify how the bench errors change the recon¬ 
structed LED spectral intensities, S (A, T). To do this, we com¬ 
pute the derivatives of the S (A, T) models with respect to all the 


identified systematics. This is done by shifting each term of the 
systematics vector in turn and re-determining the model param¬ 
eters as described in section 5.1 above. They are used to propa¬ 
gate the test-bench systematics, (see section 6 and Appendix E 
for details). 

5.4. Conclusion 

At this stage, we have modelled the emissivity of each LED as a 
smooth function that predicts the spectral intensity of the beam 
in any given direction and at any temperature. The bench statis¬ 
tical uncertainties were encoded into the covariance matrices of 
the 0 parameters. The main systematics were quantified. They 
come primarily from our calibration of the monochromator and 
from the error budget reported at NIST. How they will be propa¬ 
gated in the final analysis is discussed in section 6 and Appendix 
E. 

We now have all we need to measure the transmission of an 
imager from series of calibration frames taken with these light 
sources. This is the subject of the next section. 

6. Discussion: caiibrating broadband observations 

We have built a very stable (0.01%) light source and charac¬ 
terised it on a test bench with an accuracy of ~ 0.3 nm in wave¬ 
length and ~ 0.1% in flux. Now the question is how we use it 
to calibrate a real imager. And first of all: what do we need to 
measure in order to calibrate an imager? 

6.1. Transmissions 

Passbands are known long before the first star light hits the fo¬ 
cal plane. The transmissions of all optical components and the 
quantum efficiency curves of the detectors are measured before 
assembly. These measurements are combined to build a synthetic 
passband model; 

T(A)^gXjlX Rmirror(T) X roptics(d) X 7’fiiter(/l) X s(A) (13) 

where s(/l) is the quantum efficiency of the CCD, g is the gain 
of its readout chain, is the area of the mirror, and the other 
terms are the various transmissions and reflectivities of the op¬ 
tical elements. Note that T(A) is a dimensioned quantity; here it 
has units of ADU/y/m^. 

The absolute normalisation of T{A) varies with time and 
must be monitored: the gain g of the readout electronics may 
fluctuate by a few per-mil on timescales of a few hours; also, 
alterations of the optical surfaces (dust deposits, ageing of the 
coatings, etc.) slowly degrade the transmission of the instru¬ 
ment by as much as 5-10% per year. The resulting attenuation 
of the telescope transmission is slightly wavelength dependent. 
Therefore, the relative normalisation of the passbands with re¬ 
spect to each other may itself vary by a few percentage points 
per year. The main purpose of calibration is therefore to monitor 
these variations. In most applications, what we need to measure 
is the evolution of the relative normalisation of the passbands 
with respect to each other. 

The shape of the passbands is not expected to vary very sig¬ 
nificantly over time. This is a design requirement. However, sev¬ 
eral studies have reported evidence of slow evolution of the pass- 
band shape. For example, Doi et al. (2010) report a 30% decrease 
in the short-wavelength side of the SDSS 2.5-m w-band chan¬ 
nels, probably due to ageing of the CCD anti-reflective coat¬ 
ings. Another example can be found in Betoule et al. (2013), 
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SkyMapper g & r (distorted) 



Fig. 18. Distortion of the red and blue fronts of the SkyMapper 
g and r passbands, respectively. 


who show hints that the red fronts of the r-band and /-bands that 
equips MegaCam are about 8 -nm off with respect to the scans 
provided by the manufacturer (resulting in a disagreement of 4 
nm on the filter mean wavelength). In this case, there is indi¬ 
rect evidence that this evolution of the passbands took place in 
an early phase of the life of the instrument. There is therefore a 
strong incentive for future surveys to monitor the shape of their 
effective passbands. 


a single vector, noted 

^^blue|u 

^dj-edlu 

: 

N, 

^^bluel^ 


( 15 ) 


To constrain these unknowns, we build a model that predicts the 
flux registered on the focal plane. These predictions are built 
from the LED spectral intensity estimates S i(A, T), derived in 
section 5, and the re-parametrised passbands that we denote as 
The model simply writes as 

<Pbi = dw JS,(A, T) Tb{A,-»,) dA (16) 

where dot is the solid angle covered by a focal plane pixel or 
superpixel (about 10 “'^ sr for a single pixel, and 10 “^ sr for a 
typical 128 x 128 superpixel). Comparing measurements of the 
calibration flux performed with the imager with the predicted 
flux tphi yields constraints on 

The details of the method are described in Appendix E. 
Classically, it consists in minimising a function built from 
the model above and the flux measurements. We are careful to 
propagate the effect of the identified uncertainties to the final es¬ 
timates. The statistical uncertainties affecting the LED spectral 
intensity models are accounted for. The systematics, those re¬ 
lated to the test bench measurements, as well as those affecting 
the flux measurements on site, are incorporated into the fit as 
nuisance parameters and marginalised over. 


6.2. Constraining passbands 

Instead of measuring the filter transmissions again with a 
monochromatic beam, the DICE strategy consists in taking the 
(known) synthetic passbands as a starting point and in constrain¬ 
ing small alterations to these passbands from series of measure¬ 
ments. 

The broadband flux of an astrophysical object is primarily 
sensitive to ( 1 ) the normalisation of the passband and ( 2 ) the 
position (in wavelength) of the blue and red filter cutoffs. These 
three quantities are therefore what should be monitored in the 
long run. To do this, we alter the synthetic passband model T (A) 
presented in equation 13 to allow for a different normalisation 
and for small (potential) variations in the filter cutoffs. The latter 
is done by composing r(A) with a linear function that shifts and 
stretches (or dilates) A around the filter mean wavelength: 


A A'^ a(A - A) + j8. (14) 


In practice, we reparametrise the function above, so that it de¬ 
pends directly on the filter cutoff displacements, denoted i5/lbiue 
and ^/Ired hereafter. This way, one can shift each front easily, and 
they are essentially independent of the other, as shown in figure 
18. 

By doing so, we reduce our problem to fitting only three pa¬ 
rameters per filter: a normalisation, TV, and two filter front dis¬ 
placements. Eor an imager equipped with N (typically 5) filters, 
we end up with 3 xN calibration parameters, which we group in 


6.3. Lightweight caiibration runs 

The system has been designed to be run in routine mode every 
day or so. During a typical calibration run, each passband is sam¬ 
pled with the corresponding LEDs. Figure 3 shows that the typ¬ 
ical u{v)gnz passbands are covered by three to nine LEDs, de¬ 
pending on their extension in wavelength, with the filter cut-ons 
and cut-offs sampled by one or two LEDs"^. In Table 2, we esti¬ 
mate the number of exposures that have to be taken during typ¬ 
ical calibration runs for various designs and the time requested 
to complete these runs. Conservatively taking about one minute 
per exposure, plus 20 minutes for the overheads (filter changes, 
alignment of the source, and the telescope, etc.), we find that a 
typical run requires a little less than one hour of daytime. This 
calibration program is short enough to be run every day or so 
just after and/or before telescope operations. 

Significant gains in precision can be obtained by combin¬ 
ing several calibration runs. Indeed, if the normalisation of the 
passbands does change slightly from one night to the next, 
the position of the filter cut-offs is not expected to move in 
a measurable way overnight. As a consequence, we will typi¬ 
cally fit for one normalisation per filter and per run, but only 
one position parameter for each filter cut-off. This represents 
~ A^bands X A^runs + 2 X Abands (e.g. 60 parameters for 10 MegaCam 
calibration runs and 135 parameters for 25 runs). As shown be¬ 
low, this allows us to improve the precision on all calibration 
parameters sizeably, at the price of a slightly more complex pro¬ 
cedure. 

Except for the early version, SnDICE, which presented an under¬ 
sampled region around 700-nm, precisely at the location of the red front 
of the r-filter, and the blue front of the /-filter 
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Table 2. Typical MegaCam and SkyMapper calibration runs. 



u 

V g 

r 

i z 

Duration 
(# frames / mn) 

MegaCam 

(SnDICE) 

4 

- 7 

5 

7 4 

27/47 

SkyMapper 

(SkyDICE) 

4 

4 8 

6 

5 6 

33/53 


6.4. Simulated datasets 

To conclude, we assess the quality of the constraints that can 
be obtained from series of calibration runs of the type described 
above. The main goal here is to show that the problem is well 
constrained and to give estimates of the uncertainty budget. 
These estimates are derived from the quantitative analysis of the 
accuracy of the source test bench calibration (section 5). These 
are lower bounds, because we do not discuss the uncertainties 
associated with the imager measurements of the calibration light 
in detail. Based on preliminary analyses of the calibration frames 
taken with MegaCam and SkyMapper, we assume that they are 
of about 0.5% (uncorrelated). 

We generate synthetic flux measurements, using 
MegaCam/SkyMapper passband models as follows. We 
first simulate realistic test bench datasets, by modelling the 
“true” LED spectral intensities Stniei4,T) with Gaussians 
of width cr ~ 15-nm, normalised to a radiant intensity of 
0.1 mW/sr, and emulating the true SnDICE and SkyDICE 
coverage. We estimate the LED spectral intensities with the 
methods described in sections 4 and 5. We use the “true” 
(Gaussian) spectra to generate the broadband fluxes and the 
reconstructed spectral intensity models S (A, T) to estimate 
the calibration parameters. The calibration parameters are 
estimated following the method described in section 6.2 and 
Appendix E (using the reconstructed LED spectral intensity 
models). The method yields the full stat+syst covariance matrix 
of the calibration parameters, and we report the diagonal 
elements as a function of the number of calibration runs. 

The quality of the constrains on is a function of the quality 
of the flux measurements performed on site, but also of how we 
manage to sample each filter. We therefore study not only what 
may be obtained with the sources already in place, but also what 
could be done with optimised designs, such as those discussed 
in Appendix A. All designs involve the same number of LEDs 
(24), so that the duration of the corresponding calibration runs is 
similar to what is reported in Table 2. 


Filter cutoffs The precision on the filter cut-off positions only 
slightly depends on the number of runs that are combined into 
the calibration fit. In five to ten runs, depending on the quality 
of the cut-off coverage, it hits a systematics floor of about 0.3- 
nm, where the dominant systematics here is the precision of the 
monochromator wavelength calibration. It affects all the filter 
cut-off estimates performed with a given source. 

In figure 19 we show, for SnDICE (MegaCam), SkyDICE 
(SkyMapper), and their associated optimised designs, the un¬ 
certainties on the filter cut-off positions obtained by combining 
about 20 calibration runs. SkyDICE is very close to an optimal 
design, and it can constrain the SkyMapper filter cut-offs ex¬ 
tremely well. On the other hand, SnDICE suffers from its identi¬ 
fied lack of coverage around 700-nm. In particular, it yields only 
marginal constraints on the red cutoff of MegaCam r-band. 


E 

'O 

1 


10 '^ 


E 

1 


10 '^ 


Fig. 19. Upper panel: Uncertainties in nanometres on the blue 
(blue triangles) and red (red crosses) cut-offs of the MegaCam 
passbands, after 20 calibration runs. Grey points: the same, for 
an optimised design such as the one described in Appendix A. 
Lower panel: the same, for SkyMapper. 
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Fig. 20. Uncertainty on the MegaCam f-band normalisations, rel¬ 
ative to the r-band, as a function of the number of calibration 
runs. The uncertainty level actually depends on how we inter¬ 
pret the NIST error budget colour uncertainty (black line) or gray 
scale uncertainty (dashed red line) - see section 5.3 for details. 
The grey lines show the calibration uncertainties obtained with 
an optimised filter sampling (see Appendix A). 


Passband normalisations The uncertainties on the normalisa¬ 
tion of the MegaCam and SkyMapper passbands (relative to r) 
improve sizeably with the number of runs, as shown in figure 20. 
In one run, we are able to reach a sub-percent accuracy. After 
about ten calibration runs, the uncertainty is divided by two, be¬ 
low 0.5% in all bands. 

The level of the systematics floor depends on how we inter¬ 
pret the NIST uncertainties. In the “best case scenario” discussed 
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Fig. 21. Upper panel; uncertainty on the MegaCam passband 
normalisation (relative to r) obtained when combining 20 
SnDICE runs. Lower panel: the same, for SkyMapper. 


in section 5.3.2, in which the NIST uncertainties are interpreted 
as a grey-scale uncertainty, the impact of the NIST error budget 
is small, and we are able to reach uncertainties of 0.3% or be¬ 
low. In the other scenario, where the NIST errors are negatively 
correlated, the error floor can be significantly higher. 

Figure 21 summarises the normalisation uncertainties for 
SnDICE (MegaCam) and SkyDICE. The true uncertainties are 
somewhere in the grey band, probably closer to the “best-case 
scenario” line. This shows that neither our design choices nor 
the finite precision of the test bench measurements may limit the 
precision of the calibration that we will obtain from such mea¬ 
surements. All essentially depends on how accurately we can 
measure the calibration flux delivered on the focal plane. 

7. Conclusion 

We have presented the concept and performances of a stable 
LED-based calibrated light source, designed for the photomet¬ 
ric calibration of wide held imagers. The design has been delib¬ 
erately kept as simple as possible. The light is generated by a 
series of 24 narrow spectrum LEDs chosen to cover the imager 
spectral range. Each LED implements a point source, generat¬ 
ing a conical beam, which yields a uniform illumination of the 
focal plane. No intermediate optical surface is present between 
the LED emission zone and the telescope primary mirror. This 
makes the hardware very stable and easy to maintain in the long 
run. 

In this paper, we have described the first step of our calibra¬ 
tion program in detail (labelled O on figure 1); the transfer of the 
NIST flux scale to the light source itself. The emissivity of each 
LED is characterised by its spectral intensity, i.e. the emitted 
power per unit solid angle and per unit wavelength. The LED 
spectral intensities are measured on a spectrophotometric test 
bench, using a monochromator and a Hamamatsu S2281 photo¬ 
diode calibrated at NIST. The test bench data are combined into 


a spectral intensity model that captures the variations in the LED 
emissivity as a function of temperature (the so-called “brighter- 
cooler” and “bluer-cooler” effects). This model is then used on 
site to predict the light actually delivered by the source at a given 
temperature. 

The source is intrinsically stable at the level of < 1%. It is 
equipped with several monitoring systems to control its stabil¬ 
ity during operations. We have shown that if we monitor either 
the temperature of the LEDs or directly the light emitted by the 
LED, we can control the relative variations in the calibration 
light with a precision that ranges from a few 10 for the best 
channels to about 10“^. This means that we have built a source 
that is as stable as the silicon photodiodes used to disseminate 
the NIST flux scale. This source can then be used to calibrate 
other light detectors (e.g. astronomical imagers). 

In the last section of this introductory paper, we sketched 
the second step in our calibration program (labelled © in fig¬ 
ure 1) in order to give a complete picture of the system and its 
main use case. We described how we extract simultaneously the 
relative normalisation of the telescope passbands and the posi¬ 
tion of the passband cutoffs from a relatively small number of 
non-monochromatic calibration exposures (4 to 9, depending on 
the passband width). All the uncertainties, either those related to 
the test bench measurements or those affecting the flux measure¬ 
ments performed on site, can be propagated to the final result 
(including the non-diagonal terms). 

As an illustration, we quantified the impact of the test bench 
systematics on the calibration parameters. Regarding the Al¬ 
ter normalisations (relative to the r-band), we showed that the 
systematics floor is well below 0.5%. The exact value actually 
depends on how we interpret the NIST error budget. This is 
an important point that will have to be clarified in the future. 
Regarding the passband cutoffs, we have shown that, although 
the source is not monochromatic, the constraints we obtain are 
essentially limited by the accuracy on the wavelength calibration 
of our monochromator. 

As of today, more than 30 calibrations runs have been accu¬ 
mulated with MegaCam, about ten of them, in conjunction with 
direct observations of CALSPEC standards. The analysis of this 
dataset is ongoing. It constitutes the next step in our calibration 
programme, and it will be described in upcoming papers, cur¬ 
rently in preparation. 

Paper 2 will be devoted to analysing the calibration expo¬ 
sures, in particular determining the passband normalisations and 
cutoffs, as described in section 6 (Step © of Fig. 1). A specific 
analysis pipeline is being developed for this purpose. The pre¬ 
processing stage (e.g. bias subtractions, gain equalization, etc.) 
is extremely similar to what is implemented in a classical im¬ 
age processing framework. The downstream operations are more 
specific; the goal is to measure the irradiance delivered on the 
focal plane by an extended beam, and one has to deal with non¬ 
standard foreground effects, such as (1) the diffraction patterns 
caused by dust and optical defects in the light path, (2) the angu¬ 
lar distribution of the beam on the MegaCam interference filters, 
and (3) stray light (or ghosts) coming from reflections within the 
telescope optics. Our experience with the MegaCam imager has 
shown that ghost contamination can be as high as 20% in some 
regions of the focal plan. This effect is chromatic and can there¬ 
fore affect the relative calibration of bands. How we deal with it 
has been the central topic of this paper. 

Paper 3 will describe Step © of our calibration programme, 
i.e. the comparison between the NIST and CALSPEC flux 
scales, using the observations of the CALSPEC objects, taken in 
conjunction with the SnDICE2 calibration runs. There are two 
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main difficulties: one is the control of the atmospheric trans¬ 
mission at the time of the observations. We deal with this by 
carefully scheduling the observations, along with using a semi- 
empirical model of the atmosphere, trained on observations per¬ 
formed at the UH 2.2-m with the SNIPS instrument. (Some of 
these observations were taken in conjunction with the MegaCam 
observations of the CALSPEC standards.) The other difficulty 
resides in assessing the systematic differences between calibra¬ 
tions obtained from flat held exposures, on one hand, and point 
source measurements, on the other. Low-angle scattering of light 
within the telescope optics is the main effect that has to be un¬ 
derstood. 

For now, we have fulhlled the hrst step of our programme, 
which is to demonstrate that it is possible to build LED-based 
sources, to characterise them spectrophotometrically with an ac¬ 
curacy of a few 10 in Hux and 3 A in wavelength, and to control 
their behaviour in the long run. The DICE sources can be used 
to disseminate a NIST-based calibration to a telescope. 
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Fig.A.l. SnDICE LEDs mounted on their front-end board. 
The characteristics of all LEDs are summarised in Table A.l. 
The central LED in slot DO is used to generate the planet 
beam. The flat-top LEDs mounted in a white packaging are the 
high-intensity OSRAM Golden Dragon®. The remaining slots 
host lower intensity LEDs purchased from Marubeni and Seoul 
Semiconductors. They cover the redder wavelength range. 
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Appendix A: LED selection 

SnDICE LEDs We now give additional information about the 
calibration LEDs used for SnDICE. In the wavelength range 
400 nm < T < 630 nm, we use a series of high-intensity LEDs, 
manufactured by OSRAM Semiconductors GmbH and marketed 
under the name Golden Dragon®. The Golden Dragon® LEDs 
are flat top InGaN (up to ~ 550 nm) and InGalP (beyond 
550 ~ nm) emitters, able to deliver up to a few hundred mil¬ 
liwatts at maximum current. This is actually far more than what 
is needed. At the telescope, they are therefore operated on a low 
regime and emit about ~ 0.5 mW/sr. 

At longer wavelengths, we use AlGaAs LEDs, manufactured 
by Marubeni America Corporation (SMC series). These compo¬ 
nents are much less powerful, delivering a maximum of ten mil- 
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liwatts. They nevertheless suit our needs perfectly, as we only 
need a fraction of this illumination. 

In the UV, the LEDs that were available on the market in 
2007 were unfortunately much less powerful. We selected a 
series of components manufactured by Seoul Semiconductors. 
They all emit a maximum of a few 10“^ mW per steradian. As 
a consequence, SnDICE requires very long exposures (hundreds 
of seconds) to accumulate enough flux in the MegaCam w-band. 

Eigure A.l shows the SnDICE LED mother board. The 
Golden Dragon LEDs are the flat top components mounted in 
white packaging. The Marubeni emitters are the flat top LEDs, 
on the left. Whenever possible, we tried to choose flat top de¬ 
vices, which are Lambertian emitters with a very good approxi¬ 
mation. All the blue LEDs procured from Seoul Semiconductors 
have either a flat (D8, D21) or hemispheric (D4, D14, D16, D24) 
protection. They are nevertheless quasi-Lambertian emitters in 
the angular range we are considering (+1°). 

The wavelength coverage of the MegaCam passbands is 
shown in figure 3(a) and discussed earlier in this paper. It is 
not very satisfactory, mainly because of the gap around 700-nm, 
which leaves the red cutoff of the r-band and the blue cutoff the 
the /-band unconstrained. SnDICE was designed in 2007, and 
the LED diversity was not then what it is today. However, the 
SnDICE coverage was considered perfectly adequate for tests 
and feasibility studies. 

The types and properties of the LEDs that equip SnDICE are 
summarised in Table A.l. We report the coefficients that char¬ 
acterise the brighter-cooler and bluer-cooler relations. As can be 
seen, the mean wavelength of the LED spectra vary typically 
by 1A/°C (note however, that some LEDs seem to display an 
inverse behaviour). The importance of the flux variations as a 
function of temperature varies from LED to LED, from a frac¬ 
tion of one percentage point per degree up to 3%/°C for the UV 
LEDs. Again, one LED seems to display a different behaviour, 
at least at the nominal current at which it is operated. These re¬ 
lations are (1) well measured and (2) very reproducible, which 
is all that matters for our application. 

SkyDICE LEDs The SkyDICE source was designed about three 
years later. At that time, the LED market had literally exploded, 
and a much higher diversity of narrow spectrum LEDs was 
then available. Eigure 3(b) shows the resulting sampling of the 
SkyMapper Alters. As discussed earlier in this paper, this cover¬ 
age is much better, although the LED density in the blue part of 
the spectrum could be improved, in order to obtain tighter con¬ 
strains on the narrow u- and v-bands. 

The types and properties of the SkyDICE LEDs are sum¬ 
marised in Table A.2. Again, in the visible, we use InGaN and 
InGaAlP Golden Dragon® components. In the IR, we use a new 
family of high-intensity GaAlAs LEDs distributed by Roithner 
Lasertechnik GmbH. Einally, in the UV, down to ~ 360-nm, the 
situation has improved considerably in terms of emitted power. 
In this region, we rely on a series of InGaN LEDs from Roithner 
Lasertechnik. 

These new-generation LEDs are more powerful than their 
SnDICE counterparts, and they have to be operated at even lower 
currents. Again, the parameters of the brighter-cooler and bluer- 
cooler relations are reported in Table A.2. They are similar in 
magnitude to what was measured with SnDICE. As mentioned 
earlier in this paper, the amplitude of these effects do strongly 
depend on the LED technology and mean wavelength. 


Optimised designs The diversity of narrow spectrum LEDs 
continues to improve at a rapid pace (see upper right panel of fig¬ 
ure 3). A quick glance at the Roithner and OSRAM catalogues 
shows that there are about 50 LEDs covering the wavelength 
range 245 nm < T < 1100 nm, allowing us to sample the imager 
wavelength range with almost one LED every 10 to 20-nm. One 
may therefore wonder whether it is feasible to design a calibra¬ 
tion source that would sample the imager response in some kind 
of “optimal” way. 

There are many ways to define optimality. We seek to iden¬ 
tify the set of LEDs that would allow us to obtain the best con¬ 
straints on the model described in section 6 after one single cal¬ 
ibration run. We proceed as follows. We first identify the LED 
positions that would yield the best constraints on the Alter fronts: 
these are the positions that maximise the derivatives of the iphi 
model as a function of the Alter cutoff displacements 6Ah and 
SAr- We match these positions as precisely as possible using the 
LEDs available on the market. In general, we are able to And a 
LED within 10 nm of the optimal position. It is possible to play 
on the diversity of LEDs of the same model to reach optimality 
even better. Einally, we cover the gaps between the Alter fronts 
as densely as possible, with about one LED every 25-nm. 

Eor MegaCam, this procedure yields the coverage shown in 
the upper panel (left) of figure 3. The dotted lines show our “op¬ 
timal” positions. The blue bands show the peak emission of the 
spectra (+5% of the peak position). The LEDs available in the 
OSRAM and Roithner catalogues are indicated with arrows. In 
practice, we have seen on simulations that the optimal design 
(dotted lines) and its actual implementation (arrows) are nearly 
indistinguishable in terms of performances. Eor SkyMapper, we 
have not been able to identify a design that would give con¬ 
straints significantly better than those obtained from the current 
source. 

In section 6, we systematically compare the performances of 
the existing DICE sources with these optimised designs. In par¬ 
ticular, we show that, although SnDICE is far from optimality, 
SkyDICE is nearly optimal, except for the narrow SkyMapper 
v-band. 


Appendix B: Glitch observed on a APG2C1-760 

In section 4, we mentioned a 2% glitch observed on a Roithner- 
Lasertechnik APG2C1-760 LED during a long-duration test. In 
figure B.l we show the timelines registered by the source and 
the bench sensors. We see in the upper left-hand panel that after 
about two weeks of uninterrupted operation, the LED flux sud¬ 
denly increased by almost two percent. This increase was de¬ 
tected by the off-axis control photodiode that monitors the LED 
flux (upper right panel). However, for a reason that is not fully 
understood, it registered only a fraction of the variation (~ 1%). 

The other panels of figure B.l display the timelines of the 
source parameters, monitored along with the LED flux: the LED 
temperature, the temperature of the backend board, the LED cur¬ 
rent, and the reference tension of the electronics. No correspond¬ 
ing glitch was registered on any of these quantities. 

As of today, we have no clear explanation for this phe¬ 
nomenon. It may be due to a sudden variation in the structure 
of the LED junction. We have not been able to observe any sim¬ 
ilar event in any of our subsequent long duration tests. In any 
case, should such an event occur during operations, it would 
be detectable by the off-axis control photodiode, and the corre¬ 
sponding LED would not be used in subsequent operations (until 
recalibration of the source). 
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Table A.l. Summary of the SnDICE LED characteristics. 


Model 

Type 

^LED / hnax 

(mA) 

<4> 

(@ 25 °C) 
(nm) 

dAldT 

(nm/°C) 

Radiant intensity 
(@ 25 °C) 

(mW / sr) 

40/(0)47) 

(%/°C) 

channel 

S8D31C® 

ALGaN 

15/20 

312.5 

-0.035(027) 

0.02181(00021) 

-2.2129(0053) 

D4 

S8D34C° 

ALGaN 

15/20 

342.3 

-0.027(008) 

0.01885(00003) 

-2.1448(0101) 

D24 

T9F34C° 

ALGaN 

15/20 

387.5 

-t0.059(005) 

0.01072(00002) 

-3.0912(0166) 

D21 

T9F31C® 

ALGaN 

15/20 

388.0 

-t0.086(005) 

0.01539(00001) 

-2.1406(0050) 

D8 

S8D40® 

InGaN (?) 

17.5/350 

408.0 

-t0.029(001) 

0.49891(00018) 

-0.4522(0033) 

D14 

S8D42® 

InGaN (?) 

17.5/350 

422.2 

-t0.015(002) 

0.20007(00019) 

+0.1983(0025) 

D16 

LBWSSG* 

InGaN 

25 / 500 

466.4 

-t0.017(001) 

0.23207(00012) 

-0.2156(0005) 

DIO 

LBW5SG+ 

InGaN 

25 / 500 

468.7 

-tO.OlO(OOl) 

0.33045(00016) 

-0.2153(0016) 

D23 

LBWSSM* 

ThinGaN 

25 / 500 

474.3 

-tO.018(000) 

0.65071(00055) 

-0.2241(0007) 

D3 

LTW5SG* 

InGaN 

25 / 500 

533.5 

-0.012(000) 

0.12352(00003) 

-0.1301(0005) 

D13 

LTWSSM* 

ThinGaN 

25 / 500 

540.6 

-t0.014(001) 

0.31282(00023) 

-0.5900(0013) 

D7 

LYW5SM* 

Thinfilm InGaAlP 

25 / 500 

590.6 

+0.105(003) 

0.05893(00006) 

-0.3915(0022) 

D20 

LAW5SM* Thinfilm UnGaAlP 

25 / 500 

620.2 

+0.121(001) 

- 

- 

D19 

LRW5SM+ 

Thinfilm InGaAlP 

25 / 500 

630.3 

+0.114(001) 

0.38592(00029) 

-0.2303(0005) 

D2 

LRW5SM* 

Thinfilm InGaAlP 

25 / 500 

630.6 

+0.117(001) 

0.36542(00017) 

-0.3644(0021) 

D18 

SMC750* 

AlGaAs 

30/50 

729.4 

+0.170(002) 

0.23181(00018) 

-0.8754(0013) 

D22 

SMC735* 

AlGaAs 

30/50 

731.8 

+0.172(003) 

0.23947(00028) 

-0.4147(0008) 

D12 

SMC735*t 

AlGaAs 

25 / 500 

733.1 

+0.168(003) 

- 

- 

D6^ 

SMC750* 

AlGaAs 

30/50 

746.8 

+0.163(002) 

0.29190(00003) 

-0.2369(0009) 

D9 

SMC810*+ 

AlGaAs 

30/50 

798.4 

+0.192(002) 

- 

- 

D1 

SMC810*t 

AlGaAs 

30/50 

- 

- 

- 

- 

DIE 

SFH 4230* 

GaAs (?) 

25 / 1000 

843.3 

+0.223(001) 

0.37082(00037) 

-0.1613(0008) 

D15 

SFH 4203“* 

GaAs 

75 / 100 

942.4 

+0.308(003) 

0.69169(00091) 

-0.2485(0050) 

D17 

SFH 4203” 

GaAs 

75 / 100 

941.4 

+0.282(003) 

0.74315(00118) 

-0.5340(0020) 

D5 


O Seoul Semiconductor Co., Ltd - http: //www. acriche. com 

j: Golden Dragon® series, OSRAM Opto Semiconductors GmbH - http://www. osram-os. com 
★ Marubeni America Corporation - http: //tech- led. com 
X OSRAM Opto Semiconductors GmbH - http://www. osram-os. com 
t dead channel 


Appendix C: Spectroscopic caiibration bench 

Test bench setup and operations The test bench setup is pre¬ 
sented in figure 8. The monochromator is operated in a slightly 
non-standard way, in the sense that the LEDs are not placed just 
in front of the entrance slit, but at a distance of 270 mm. Since 
the LED emissive zone is very small, the light that enters the 
monochromator is a quasi-pencil beam, with an angle distribu¬ 
tion less than 0.2°wide. 

The value of the mean angle of the incident beam is defined 
by the relative positioning of the LED and the entrance slit. We 
(conservatively) estimate a positioning error of the LED head 
with respect to the monochromator of about 1 mm. This trans¬ 
lates into an uncertainty on the incident beam angle of 0.2°. 

The output is a larger beam, convolved with the monochro¬ 
mator response. At the level of the photodiode, it produces an 
illumination of width 4 mm and height 6 mm with an extended 
penumbra zone. It is almost entirely contained in the S2281 
NIST photodiode used to measure it. Since we were concerned 
that there could be chromatic flux losses affecting the measure¬ 
ments, the beam was scanned spatially and in wavelength in or¬ 
der to check for diffuse light and parasitic reflexions. We found 
that the flux loss is essentially constant (0.5%) as a function of 
wavelength in the 400 nm < 4 < 800 nm range. It then increases 
up to 1% for redder wavelengths. The LED SED being narrow, 
the shapes of the measured spectra are not affected by this effect. 

The size of the entrance and output slits define the 
monochromator wavelength passband -through the dispersion 
relation. The latter has been remeasured by studying images 


of monochromatic sources with various slits apertures. We find 
3.2 ± 0.11 nm per millimetre aperture in accordance with the 
manufacturer indications. We set the slit openings to 625 pm and 
the monochromator step to 2 nm for optimal sampling, during 
the full calibration process. This choice is a trade-off between 
fineness of sampling and S/N (which is critical for UV LEDs). 

Apart from the grating and picoammeter range, which dif¬ 
fer for each LED, the measurement protocol is the same for all 
calibration channels. The temperature is set to the desired value, 
and we let the bench thermalise. Each LED is positioned in front 
of the entrance slit, and a fast wavelength scan is performed in 
order to check the SED extension. Each spectrum is then sam¬ 
pled every 2 nm. At each step, we alternate measurements of the 
dark current and the LED flux. Spectra are obtained at least at 
six different temperatures ranging from 0 °C to 25 °C. 


Wavelength Calibration The monochromator does not come 
with a temperature-dependent wavelength calibration. During 
the tests, the monochromator temperature has been found to 
vary between 5°C and 25°C. (The bench is never fully ther- 
malised, and the monochromator is always slightly warmer 
than the LEDs.) It is therefore essential to check for potential 
temperature-dependent effects. The absolute wavelength calibra¬ 
tion of the monochromator was determined at 25°C, 18°C and 
13°C using two calibration lamps: one Cd lamp and one Hg 
lamp. We have found a small wavelength offset, AT, that varies 
slightly as a function of the temperature of the monochromator. 
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Table A.2. Summary of the SkyDICE LED characteristics. 


Model 

Type 

^LEd! hna.x 

(mA) 

<4) 

(@ 25 °C) 
(nm) 

dAldT 

(nm/°C) 

Radiant intensity 
(@ 25 °C) 
(mW / sr) 

d9l{(t>dT) 

(%/°C) 

channel 

UVTOP315-FW-T039* 

AlGaN 

12.2/20 

320.7 

-t0.001(014) 

0.0866(0002) 

-0.330(014) 

D8 

UVTOP335-FW-T039* 

AlGaN 

12.2120 

340.4 

-t0.044(025) 

0.1081(0002) 

-0.665(009) 

D21 

APG2C1-365-S * 

InGaN 

21.4/350 

368.3 

-t0.027(006) 

0.6625(0006) 

-1.534(004) 

D2 

APG2C1-385*' 

InGaN 

-/350 

- 

- 

- 

- 

D23t 

APG2C1-395* 

InGaN 

10.7/350 

396.7 

-t0.020(019) 

0.9606(0006) 

-1.003(003) 

D3 

APG2C1-420* 

InGaN 

2.1/350 

417.0 

-t0.015(016) 

0.2486(0003) 

-0.535(005) 

D22 

LD WSAM* 

Thin GaN 

3.1/500 

452.6 

-t0.013(009) 

1.1194(0010) 

-0.117(005) 

D4 

LB W5SM" 

Thin GaN? 

3.1/500 

466.6 

-t0.016(019) 

1.6429(0027) 

-0.080(009) 

D17 

LV W5AM+ 

Thin GaN? 

3.1/500 

515.9 

-0.004(012) 

0.6586(0001) 

+0.256(001) 

D5 

LT W5SMt 

Thin GaN 

3.1/500 

546.8 

-t0.010(008) 

0.6469(0002) 

+0.016(001) 

D24 

LT W5AMt 

Thin GaN 

3.1/500 

528.8 

-0.025(020) 

0.3338(0001) 

-0.372(001) 

D6 

LY W5SM* 

Thin InGaAlP 

15.3/500 

590.8 

-tO. 102(009) 

0.3544(0004) 

-0.268(007) 

D20 

LA W5SMt 

Thin InGaAlP 

3.1/500 

622.6 

-tO. 104(020) 

0.2292(0001) 

-0.515(003) 

D7 

APG2C1-660* 

GaAlAs 

15.3/500 

656.8 

■+0.144(009) 

0.5411(0004) 

-0.418(004) 

D18 

APG2C1-690* 

GaAlAs 

3.1/500 

687.6 

+0.155(020) 

0.1697(0001) 

-0.433(002) 

D1 

APG2C1-720* 

GaAlAs 

3.1/500 

716.5 

+0.153(020) 

0.2319(0001) 

-0.434(001) 

D16 

APG2C1-760* 

GaAlAs 

3.1/500 

759.9 

+0.176(017) 

0.2686(0001) 

-0.166(001) 

D15 

APG2C1-810* 

GaAlAs 

3.1/500 

806.4 

+0.177(031) 

0.2785(0001) 

-0.197(002) 

DIO 

APG2C1-830* 

GaAlAs 

3.1/500 

828.2 

+0.203(011) 

0.3401(0003) 

-0.352(004) 

D19 

APG2C1-850* 

GaAlAs 

3.1/500 

845.7 

+0.189(025) 

0.2519(0002) 

-0.297(004) 

D12 

SFH42E 

Thin InGaAlP 

6.1/100 

730.8 

+0.153(009) 

0.4493(0001) 

-0.566(001) 

D13 

APG2C1-905* 

GaAlAs 

3.1/500 

913.6 

+0.133(087) 

0.1497(0001) 

-0.437(004) 

D9 

APG2C1-940* 

GaAlAs 

15.3/500 

950.5 

+0.142(016) 

0.5216(0007) 

-0.485(007) 

D14 

APG2C1-970* 

GaAlAs 

3.1/500 

951.4 

+0.272(059) 

0.1441(0001) 

-0.343(003) 

Dll 


* Roithner Lasertechnik - http: //www. roithner- laser. com/. 

j: Golden Dragon® series, OSRAM Opto Semiconductors GmbH - http://www. osram-os. com/osram_os/en/. 
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Table C.l. Absolute wavelength calibration offsets. 


Grating 

# 

AAgranng @ Tl=25°C 
(nm) 

Pgrating ~ dKXj(IT 

(°C) 

1 

+0.154 + 0.246 

-0.02499 ± 0.0023 

2 

+0.020 ± 0.059 

-0.02526 + 0.0016 

3 

-0.008 ± 0.096 

-0.02654 + 0.0022 


Tm- 

A A = AAgrating + Pgrating X (Tm - . (C.l) 

The values of AAgrating Pgrating are reported in Table C.l. 

This correction was applied to all wavelength measurements. 
The final wavelength calibration error is classically derived by 
propagating the uncertainties on AAgratings and Pgrating and on the 
monochromator temperature (ctt^, ~ 1°C). We find that it never 
exceeds lA (for UV LEDs). 

Transmission The monochromator transmission is a rapidly 
varying function of A. It cannot be considered to be flat on the 
scale of one single LED spectrum, so it has to be measured pre¬ 
cisely. To perform this measurement, we placed a 800 jim wide 
mask in front of the monochromator entrance slit, at a distance 
of about 1 cm. The latter was opened slightly wider than in nor¬ 
mal operations (1mm instead of 625 pm), so that it contains a 
large, constant fraction of the beam. 

The transmission of the monochromator at a given wave¬ 
length A can then be estimated using a LED emitting in that 


wavelength range, by measuring the ratio between the flux mea¬ 
sured with the monochromator <p(A, LED) and the flux without 


it; 


T(A)^ 


<t>(A,LED) ^ T{A) rj(A) Sled(A) 
4^tot{AED) J" ti(A) S i^££)(A)dA 


(C.2) 


The monochromator transmission can be modelled knowing the 
groove spacing (1200 g/mm). It is a simple function of the blaze 
angle (e), the Ebert angle (a), and the focal distance /; 


T(A) oc sine 


7rcos(e) 




2 cos 




sin(e). 


A2 


4/2 COs2(|) 


(C.3) 

Eigure C.l shows the measurements performed over the full 
wavelength range of grating #2 (500nm). It shows excellent 
agreement of the model and the observations. Eor example, for 
this grating, a fit of the model yields Otiaze - 18.1 + 0.28° and 
ffehert - 16.6 + 1°, in remarkable agreement with the constructor 
data. 

Our transmission model does not include the reflectivities of 
the four mirrors that are inside the monochromator. Inspecting 
the residuals of figure C.l, we see a dip around 550 - 600 nm 
that could be the signature of Al reflectivity. Incorporating an 
Al reflectivity function (at the fourth power because there are 
four mirrors) into the transmission model does not improve the 
fit significantly. We find, however, that correcting locally for the 
dip has a negligible impact on the spectrum shape (less than lA 
for the mean wavelength). We therefore use the model as it is 
and consider that the uncertainties it carries essentially have no 
impact on the spectrum shapes. 
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Fig.B.l. Time series registered during a long duration test of a a Roithner-Lasertechnik APG2C1-760 LED mounted on a DICE 
source. After 18 days of operations, the flux registered by the NIST photodiode (upper left panel) suddenly increased by almost 
2%. Only a fraction of this variation was detected by the internal off-axis control photodiode (upper right panel). The other panels 
display the timeline of the other quantities monitored by the backend electronics: the LED temperature (middle left panel), the 
backend board temperature (middle right panel), the LED current (lower right panel), and the board reference tension (lower right). 
No corresponding glitch can be observed on these quantities. 



A (nm) 

Fig. C.l. Top panel: Model (dark line) and measurements (trian¬ 
gles) of monochromator grating #2. Bottom panel: fit residuals 
(same units as top panel). 


Appendix D: Isotropy of the LED illumination 

In this section, we report on a study conducted to look for a 
possible dependence between the LED spectra and the direction 
of emission, u. We know already that the beam intensities vary 


by about 1-2% as a function of u. These non-uniformities are 
measured, as described in sections 3.4 and 4.3, and the resulting 
beam maps are included in the source emissivity model. Now, 
we want to verify that the shape of the source spectrum also does 
not depend on the direction of emission. 

In principle, we do not expect a large effect. Indeed, the pro¬ 
cess that generates the LED light takes place everywhere in the 
junction. Eor small angles, the photon generation and extraction 
mechanisms are the same. However, we cannot exclude effects 
induced by the LED packaging or possibly interferences taking 
place within the junction. 

This work was initiated after the DICE sources had been in¬ 
stalled on site. As a consequence, it was not performed on the 
real sources, but on a small model built by gluing a selection of 
LEDs on a radiator, and connecting them to a laboratory current 
source. This system is slightly less stable in flux than the real 
DICE sources (0.1 % stability). However, it is stable enough for 
our purpose. 

Our selection of LEDs is listed in Table D.l. They are cho¬ 
sen so as to cover the full spectral range of DICE (3800A < T < 
9500A). Also, we tried to probe as many technologies and pack¬ 
aging types as possible. All these LEDs are narrow-spectrum 
emitters with typical smooth spectra (figure D.l). There is one 
exception: the reddest LED (SFH4203) that displays sharp fea¬ 
tures probably due to fringing within the substrate. Eor this LED, 
we expect a more pronounced dependence of the spectrum shape 
on the emission angle. 

The source model is mounted on a support that can rotate 
around a vertical axis, and the LED light is injected at various 
angles wwith respect to the direction of normal emission into 
the Digikrom DK240 monochromator that equips the bench. A 
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Table D.l. LED models tested for isotropy. 


LED model 

Type 

Packaging 

(A) 

[nm] 

Tp (AT) 

°K 

^ ’runs 

Aangle 

[“] 

^spectra 

duration 

(days) 

APG2C1-395* 

InGaN 

Hemispheric 

400.3 

301.5 (1.3) 

5 

[-7.5,7.51 

126 

4.4 

LT W5SM-1XKX-36+ 

Thin GaN 

Flat top 

536.9 

301.8(3.9) 

13 

[-8.0, 8.0] 

465 

25.1 

LA W5SM-GYHZ-35+ 

“ (?) 

Flat top 

586.2 

303.4 (5.6) 

6 

[-15., 15.] 

298 

5.7 

LY WSSM* 

Thinfilm InGaAlP 

Flat top 

591.2 

301.3 (2.8) 

7 

[-6.8,8.21 

405 

7.1 

APG2C1-810* 

GaAlAs 

Hemispheric 

807.0 

300.9 (0.7) 

5 

[-7.5,7.51 

243 

8.4 

SFH4203" 

GaAs 

Flat top 

944.3 

300.2 (3.6) 

7 

[-7.5,7.51 

1702 

75.0 


* Roithner Lasertechnik - http: //www. roithner- laser. com/. 

j: Golden Dragon® series, OSRAM Opto Semiconductors GmbH - http://www. osram-os. com 
X OSRAM Opto Semiconductors GmbH - http://www. osram-os. com 


photodiode, read out with a Keithly picoammeter, is placed at 
the exit slit of the monochromator. 

Information about the dataset accumulated with this setup 
can be found in Table D.l. For each LED, series of consecu¬ 
tive spectra are taken at four to eight different angles, spanning a 
range of about ±8 degrees, around the direction of normal emis¬ 
sion. This range is chosen to be deliberately larger than the 1° 
field of view of MegaCam in order to enhance our chances of 
detecting an effect. This dataset represents hundreds of spectra 
per LED, taken over the course of several weeks. Incidentally, it 
allowed us to check the stability of the LED spectra over long 
periods of time. 

We call run a series of consecutive spectra taken at a same 
angle. Each single run is analysed independently from the others. 
The only piece of information shared from one run to another 
is the measurement of the dark current, performed during dedi¬ 
cated data-gathering sessions. The analysis consists in extracting 
from each run an average spectrum template S (A), modelled in 
practice as a series of 1-nm-wide bins. If the temperature was 
constant with a precision of a tenth of a degree, its effect on the 
spectral shapes (the cooler-bluer effect) could be neglected, and 
the model could be obtained simply by averaging the flux mea¬ 
surements in each single bin. 

In fact, the temperature has been found to vary by up to a few 
degrees during each single run. As a consequence, spectra taken 
at different angles were taken at slightly different temperatures. 
This induces spectrum-shape variations that are independent of 
the effect we are looking for. To account for this, we fit the fol¬ 
lowing (temperature dependent) model on each run; 

i(A, T) ^ S {A)+ {T + 6T- Tp) x — (T) + Mark 

oT 

where i(A, T) is the measured photocurrent, Mark the dark cur¬ 
rent, T the bench temperature reported by the bench DAQ at 
the time of measurement, Tp a median pivot temperature com¬ 
mon to all the runs performed with a given LED and reported 
in Table D. 1, an offset introduced to account for the fact that 
T lags by a fraction of a degree on the real LED temperature. 
Here, S (T) is the average LED spectrum at the pivot temperature, 
dS IdT(A) corresponds to the (binned) derivatives of the spec¬ 
trum with respect to temperature, and /dark is measured indepen¬ 
dently. Therefore, the parameters that are fitted are 5, 55 /dT, 
and the dT offsets. This model has slight degeneracies, which 
we break by adding constraints to the of the form {dTl&f' 
with (T ~ 0.1 °C. 

In figure D. 1 , we show the reconstructed LED models at their 
respective pivot temperatures for all the angles studied. For all 
LEDs that emit in the visible, the spectra taken at different angles 


are virtually indistinguishable. In the corresponding insets, we 
show the variations in the LED mean wavelengths (A), which 
are found to be stable at better than lA over an angular range 
that is much greater than the acceptance of MegaCam. 

The situation is different for the IR emitter, (SFH4203). We 
see “spectral features” whose positions and relative amplitudes 
depend on the angle of emission. The question is how this affects 
the passband measurements. 

Following the analysis sketched in section 6, the bench spec¬ 
tra are combined with a (parametrised) passband model to pre¬ 
dict broadband fluxes, 


d) = 


f 


S(A)T(A)dA, 


(D.l) 


which are then compared with the fluxes measured directly on 
the calibration frames. An estimate of how much the LED broad¬ 
band flux varies over a 1° angle (which is precisely the size of 
the MegaCam focal plane) is 


1 50 ) 

d) da ’ 


(D.2) 


where a is the direction of emission. 

We use a rectangular filter to compute the broadband LED 
fluxes from the spectrum models fitted on each run. Here, A is 
computed numerically from the broadband fluxes at a given an¬ 
gle and the broadband flux at a reference angle. We vary the 
filter cutoff with respect to the LED spectrum in order to explore 
all the configurations. In figure D.2, we show the value of A as 
a function of the red filter cutoff for two of our LEDs. When 
the red filter cutoff moves towards the redder wavelengths, the 
LED spectrum is fully encompassed by the filter shape, and the 
broadband flux becomes independent of the LED spectral shape. 
On the other hand, when the cutoff moves towards bluer wave¬ 
lengths, we explore the region where only a fraction of the LED 
light is integrated by the filter. This is the regime that constrains 
the filter cutoff positions. 

For a normal LED (left), we see that A is of the order of 10 
Only in the regime where a small fraction of the LED spectrum 
tail overlaps the filter do we observe small (but not significant) 
deviations around 0.1%. Conversely, for the IR-LED SFH4203, 
we observe that the spectrum shape variations induce variations 
in the broadband fluxes of about 0.5%. We conclude that except 
for one LED model, for which an effect was easily detected, no 
emission-angle dependent variation in the LED spectra could be 
seen over a range of angles of the order of +8°. 
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A [nm] 

(b) LT W5SM-JXKX-36 (301.8 K) 






Fig.D.l. All panels show the superposition the LED spectrum models, htted on each run. Each model therefore corresponds to a 
specific angle of incidence. Except for the IR emitter (SEH4203), the spectra at various angles can barely be distinguished from each 
other. The insets show the spectrum average wavelength (A) - f AS (A)dAI J S (A)dA computed on each run. This quantity varies by 
less than lA(peak-to-peak) over a ~ 15° range. Except for SEH4203, these variations are likely to be due to residual temperature 
effects, which not well accounted for by the crude model used for this study. 
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-'‘cutoff [nm] 

(a) LT W5SM-JXKX-36 (301.8 K) 



Fig.D.2. Broadband flux sensitivity (A) to emission-angle-dependent spectrum variations as a function of the red-cutoff of a rectan¬ 
gular broadband Alter. 


Appendix E: Constraining passbands (detaiis) 

We now explain the estimation of the calibration parameters 
from series of DICE measurements (section 6). The predicted 
flux registered on the focal plane is given by equation 16, which 
can be rewritten in matrix form: 

E (E.l) 

where 6co is the solid angle covered by the focal plane pixel or 
superpixel, are the parameters of the LED spectral intensity 
model S/(A), and E is a vector whose components are the inte¬ 
grals of the basis functions defined in section 5, convolved with 
the telescope transmission Tt{A): 

= f T) Th(A,&,) dA. (E.2) 

The calibration parameters are determined by minimizing 
a built from the measurements of the LED calibration light 
observed through the different telescope passbands (^/,/, and from 
the corresponding predictions Cpti- If we note R the vector of 
residuals, this writes as 

- -^^(Cstat + Cled + Cmodel) ^ R, (E.3) 

where Cstat is the covariance matrix of the flux measurement 
uncertainties. It will be discussed in detail in the next paper of 
the series. Here, Cied accounts for the LED intrinsic variabilities 
and Cmodei accounts for the model uncertainties. The elements of 
Cmodei are propagated from the covariance matrix of the spectral 
intensity model parameter, Cg: cov((f)hi, (pb'i) - 

Since ipu is not a linear function of we linearise it: 

ni = 6oj 0j^(Eo + T 8&,) (E.4) 

at each minimisation step. Here, So is the “current” value of the 
2 vector at a given step, and F is a matrix containing its deriva¬ 
tives with respect to the calibration parameters. Both are com¬ 
puted numerically, the full minimisation taking a little less than 
five seconds on a laptop. 


Systematics The systematics we have to consider come from 
two main sources. Eirst, we have the uncertainties that affect 
the test-bench measurements. Those have been discussed in sec¬ 
tion 5 and summarised as seven parameters listed in Table 1. 
We group them into a vector ijb, which comes along with a (di¬ 
agonal) covariance matrix Ci,- These bench systematics have an 
impact on our estimates of the LED spectral intensities, Si{A), 
and therefore on our broadband flux predictions, (pbi- 

Second, the calibration measurements performed with the 
imager are themselves affected by several systematics. They in¬ 
clude uncertainties on (1) the LED temperature, (2) the position¬ 
ing of the source with respect to the telescope, (3) systematic 
drifts of the readout electronics during the data-taking sequence, 
and (4) contamination of the calibration frames by stray light, 
in particular the ghosts generated by parasitic reflections on the 
optical surfaces. The most problematic ones are those which 
are wavelength dependent, in particular the ghost contamina¬ 
tion, which is higher when illuminating the telescope with LEDs 
close to the Alter cutoff. These contributions will be discussed 
and quantified in detail in Paper II. We will leave them aside for 
now. Our purpose is to describe the propagation method and to 
evaluate the impact of the bench systematics. 

All the contributions listed above are included as nuisance 
parameters into the calibration fit described in section 6.2, and 
marginalised over, using their uncertainty estimates (encoded in 
matrix Ct) as priors in the x^ - this is the most direct and exact 
way to compute their contribution to the total error budget. This 
means that we minimise the following ;t'^: 

+ Cied + (E.5) 

Again, the fit is non-linear, and the model is linearised at 
each step, as follows: 

iPbi = (So + Y 68) + 6iib, (E.6) 

which is the equivalent of equation E.4, with the systematics. 
The At with systematics is very fast, and takes seconds on a stan¬ 
dard laptop. 
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a, a; s\ m, a; a, a; ja, a; a. 



Fig.E.l. Correlation matrix yielded by the passband (syst. in¬ 
cluded). We note the strong positive correlations between the 
filter fronts. This is because the all uncertainties on the filter 
front positions are dominated by the wavelength calibration of 
the bench monochromator. 


ASIC. Each module is connected with an analogue link to a 
small backend electronics whose main purpose is to digitise the 
photodiode current and store it in a buffer. Thanks to their small 
size (100mm x 100mm x 15mm), these modules can be hxed di¬ 
rectly on the telescope, close enough to the calibration beam to 
intercept it. One experimental CLAP module has been installed 
for tests inside MegaPrime in front of the hlter jukebox. Besides 
demonstrating that this miniaturised low-cost monitoring solu¬ 
tion works, our goal is to investigate whether we could gain in 
precision by polarising the photodiode in reverse and operating it 
in photoelectric mode, as recommended by all photodiode mak¬ 
ers. This would have the advantage of increasing the sensitivity 
of the detector, at the expense of having a stronger dark current. 

The S3477-04 is slightly smaller than the S2281 
(5.8 mm x 5.8 mm) and has a similar sensitivity. It is cou¬ 
pled to a small Peltier effect, which allows one to operate it 
around -20 °C, hence reducing the dark current very signifi¬ 
cantly when the photodiode is polarised in reverse. NIST does 
not provide a calibration service for this type of detector, so 
we have to calibrate it ourselves, using a NIST photodiode as a 
primary standard and a DICE apparatus as an intermediate light 
source. 


Covariances The test bench systematics are shared by all 
LEDs. As a result, we expect them to introduce sizeable off- 
diagonal terms in the covariance matrix of the calibration pa¬ 
rameters. The calibration ht described above yields the full 
(stat-Hsyst) covariance matrix of the calibration parameters, 

In figure E.l, we show a typical correlation matrix obtained from 
the ht of ten (simulated) calibration runs. The matrix elements 
labelled Nx are related to the hlter normalisation (relative to the 
r-hlter normalisation). The elements labelled 6Ax are the hlter 
cutoff displacements. The matrix presented here corresponds to 
the “best case scenario” presented in section 5.3.2 and to where 
the NIST uncertainties are all fully positively correlated (best 
case scenario). The band-to-band correlations of the hlter nor¬ 
malisations are essentially negligible, while the hlter’s front dis¬ 
placements are all positively correlated. Indeed, they all share 
the uncertainty on the monochromator wavelength calibration, 
which is the dominant contribution to their error budget. 


Appendix F: Cooled Large Area Photodiodes 

The time between two test-bench recalibrations of the illumina¬ 
tion system is expected to be long - several years at least. For 
this reason, it is recommended to install a calibrated photodi¬ 
ode on site that can measure the light actually delivered by the 
device and monitor any unexpected drift in the the calibration 
beam. A solution would be to install in the dome, or even directly 
on the telescope, a calibrated photodiode procured from NIST 
coupled with a picoammeter. The typical irradiance of the cali¬ 
bration beam close to the focal plane is about 0.5 nW cm“^. A 
typical photodiode that can be procured from NIST, such as the 
Hamamatsu S2281 with an active area of lcm“^ and a sensitiv¬ 
ity of 0.3 A at 555 nm, would generate a current of 150 pA, 

which can be easily measured using a commercial picoammeter. 
This solution is nevertheless expensive and not very practical, 
because the photodiode and picoammeter are quite sizeable. 

The DICE team has developed miniaturised modules com¬ 
prising a Hamamatsu S3477-04 photodiode coupled with an 
ultra-low noise current amplifier implemented as a custom-made 
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